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1. Introduction 

1.1. Statement of the results. In this paper we analyze simple models of classical chaotic 
open systems and of their quantizations. They provide a numerical confirmation of the 
fractal Weyl law for the density of quantum resonances of such systems. The exponent in 
that law is related to the dimension of the classical repeller of the system. In a simplified 
model, a rigorous argument gives the full resonance spectrum, which satisfies the fractal 
Weyl law. Our model is similar to models recently studied in atomic and mesoscopic 
physics (see £12.41 below ). Before stating the main result we remark that in this paper we 
use mathematicians' notation h for what the physicists call h. That is partly to stress that 
our h is a small parameter in asymptotic analysis, not necessarily interpreted as the Planck 
constant. 

Theorem 1. There exist families of symplectic relations, B C T 2 x T 2 , and of their 
(subunitary) quantization, e /2(C ), iV = {2ttK)~ 1 , such that 

#{A e Spec(Bfc) : |A| > r} = c(r) h~ v + o(h~ v ) , 
r>0, h = h k = (27iD k y 1 -> 0, k -> oo , 
i/ = dim(r_(B)nW + (S)), C (r) = (27r)^ X[0 ro(5)] (r), < r (B) < 1 , 

where the integer parameter D > 1 depends on B. The set V-(B) C T 2 is the forward 
trapped set of B and W + (B) is the unstable manifold of B at any point of T -(B). The 
eigenvalues are counted with multiplicities. 

In the model discussed in detail we took D = 3. The asymptotics are actually much 
more precise and include uniform angular distribution (see Prop. IB3|) . The resonances lie 
on a lattice, and some of this structure is also seen in numerically computed more generic 
situations (some numerical results have been presented in [101 EH1 139J)- Each symplectic 
relation B (or "multivalued symplectic map") is defined together with the probabilities, 
for any point, to be mapped to each of its images: B thus represents a certain stochastic 
process. The quantizations B^ quantize the relations together with their jump probabilities 
in the precise sense given in £ 14.41 
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In the models used in Theorem ^ we can compute the conductance and the shot noise 
power (or the closely related Fano factor) — see £ J2.4.3l and references given there for physics 
background and ^jBlfor precise definitions. 

Theorem 2. Suppose that the models in Theorem^ have the openings consisting of two 
"leads" of equal width (see $6.1 for a detailed description) , so that each lead carries the 
same number, M(h) ~ h~ l , of scattering channels. Then, the quantum conductance ()6.2|) 
between the two leads satisfies 



where the exponent v is the same as in Theorem^ 

The theorem should be interpreted as follows. In (jl.lj) we see that for a model of 
scattering through a chaotic cavity, approximately one half of the scattering channels get 
transmitted from one lead to the other, the other half being reflected back (this is natural 
and well known). Asymptotics in (jl.2j) are more interesting. We see that the fractal 
Weyl law, h~ u , appears in the expression for the Fano factor. In the interpretation of the 
Fano factor in terms of "shot noise" (see §2.4. 3j) . 11/80 gives the average "shot noise" per 
"nonclassical transmission channel". This number is close to the random matrix theory 
prediction for this quantity, namely 1/8 (213 EZj. I n fact, had (|1.2|) come from a physical 
experiment rather than an asymptotic computation, it would be regarded as being in a 
very good agreement with random matrix theory 1 . 

Much of the paper is devoted to rigorous definitions of the objects appearing in the 
statements of the two theorems. In this section we give some general indications, with 
detailed references to previous works appearing below. 

We consider the two-torus T 2 = [0, 1) x [0, 1) as our classical phase (with coordinates 
P = {QiP))- Classical observables are functions on T 2 and classical dynamics is given in 
terms of relations, B C T 2 x T 2 , which are unions of truncated graphs of symplectic (area 
and orientation preserving) maps T 2 —>■ T 2 . An example is given by the baker's relation 



This is a "rectangular horseshoe" modeling a Poincare map of a chaotic open system: some 
points (here {p : l/3<g<2/3}) are thrown out "to infinity" at each iteration. 



We are grateful to Yan Fyodorov for this amusing comment. 





The Fano factor ()6.3|) is given by 






p' = p/3, 0<g<l/3 
j/ = (p + 2)/3, 2/3<g<l. 
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For relations such as B we can define the forward and backward trapped sets (see (|2.4jl 
for the definition in the case of flows): 

per_ ^ 3{pj}£Lo , po = p, {Pi\ Pj-i) e B , j > o , 
p g r + ^ 3{p i }° = _ 00 , po = p , (pf, pj-i) eB, j < o . 

In the example (jl.3j) . L = Cx [0, 1), T + = [0, 1) x C, where C is the usual |— Cantor set. 

We also define the trapped set K = T + fl T_ and, at points of f^, the stable and unstable 
manifolds, W±. In the case of the above baker's relation, 

1 log 2 

v = dimT_ n W + = - dimK = dimT, nlf_ = — — , 
+ 2 + log 3' 

but for general (possibly multivalued) relations these equalities do not hold. 
A quantization (in the sense made rigorous in < J4.5j) of B is given by 

/ Fn/3 \ 
(1.4) B h = F* N \ , h = (2vriV)- 1 , 3\N , 

\ F N/3 J 

where Tm is the discrete Fourier transform on C A/ . 
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fc = 4 
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53 


45 


40 


33 


33 


30 
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jfe = 5 


124 
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fc = 6 


237 


196 


161 


150 


142 


131 


128 


12 



Table 1. Number of eigenvalues of Bh in the regions { |A| > r }, for 27ih = 
1/N, N given by powers of 3. 



Table 121 shows the analogies between the eigenvalues of this subunitary quantum map 
and the resonances of a Schrodinger operator for a scattering situation (see £ 12. lj) . 

For Bh given by (jl.4J) we are unable to prove the fractal Weyl law presented in the 
last line of Table |2l but numerical results strongly support its validity jJU] . A striking 
illustration is given by tripling N, in which case the number of eigenvalues approximately 
doubles, in agreement with the fractal Weyl law — see Table ^ 
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The family of subunitary quantum maps in the main theorem is obtained by simplifying 
Bh, and is described explicitly in (|5.2|) . It is a quantization of a more complicated multival- 
ued relation for which r + = T 2 , L = Cx [0, 1), and T_ n W + ~ C — see Proposition 15.11 
Theorem ^ follows from the more precise Proposition 15.41 

1.2. Organization of the paper. In §21 we present related results from recent mathe- 
matical, numerical, and physics literature. In particular, in < J2.4.3I we give the physical 
motivation for the objects appearing in Theorem |21 above. §21 is devoted to the review of 
classical dynamics used in our models, stressing the dynamics of open baker's relations. 

In 21 we fi rs t review the quantization of tori. We assume the knowledge of semiclas- 
sical quantization in T*M n (pseudodifferential operators) but otherwise the presentation 
is self contained. The definitions of Lagrangian states associated to smooth and singular 
Lagrangian submanifolds is based on the ideas of Guillemin, Hormander, Melrose, and 
Uhlmann in micro local analysis but, partly due to technical differences, we give direct 
proofs of the properties we need in this paper. These properties are used to analyze the 
quantizations of the baker's relation coming from the work of Balazs, Voros, Saraceno and 
Vallejos. 

Numerical results for the (usual) quantization of the open baker's relation have been 
presented in |40| . we briefly summarize them in ^4.61 In §S] we discuss the toy model 
Bh, with two different interpretations. That section contains the proof of Theorem ^ 
Finally in ^we give precise definitions of objects appearing in Theorem Eland in a lengthy 
computation we give its proof. 

Acknowledgments. We are grateful to Christof Thiele and Terry Tao for pointing out the 
'Walsh" interpretation of our toy model, and to the anonymous referee for his comments. 
The first author thanks Marcos Saraceno for his insights on that model, and Andre Voros 
for interesting questions. He is also grateful to UC Berkeley for the hospitality in April 
2004. Generous support of both authors by the National Science Foundation under the 
grant DMS-0200732 is also gratefully acknowledged. 

2. Motivation and background 

In this section we discuss motivating topics in mathematics and theoretical physics, and 
survey related results. 

2.1. Schrodinger operators. The original motivation comes from the study of resonances 
in potential scattering. The simplest case is given by considering the following quantum 
Hamiltonian: 

(2.1) H = -h 2 A + V(q) , V G C c °°(M n ; R) . 

By assuming that the potential vanishes near infinity and that it is infinitely differentiable, 
we eliminate the need for technical assumptions - - see [22] and for more general 
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settings, in the analytic and C°° categories respectively. For instance, as noted in (5UJ 
(c.32)-(c.33)] the theory of [22] applies to arbitrary homogeneous polynomial potentials at 
nondegenerate energy levels. 



h -> 


N = (2nh)- 1 -> 00 


X exp(-it(-h 2 A + V)/h)x, t>0, 
X a cut-off on the interaction region 


B*, t = 0, 1, - - - 
-B/j a subunitary matrix 


e -itz/h^ z a resonance of H = —h 2 A + V 


A*, A an eigenvalue of B h G £(C ) 


z e [E - h,E + h] -z[0,7/i] 


1 > |A| > r > 


#{z e[E-h,E + h]-i[0, 7/1]} ~ C( 7 ) h~" E 


#{A, |A| > r} ~ C(r) JV" 



Table 2. Analogies between Schrodinger propagators and open quantum maps. 



Before discussing open systems we recall the well known results for closed systems, ob- 
tained for instance by considering H above on a bounded domain Q C M n and imposing 
a self-adjoint boundary condition at dQ (Dirichlet or Neumann). Then the spectrum, 
Spec(H), of H is discrete and, at a non-degenerate energy level E its density is described 
by the celebrated Weyl law: 

(2.2) # { Spec(B) n [E - 6, E + 5} } = -\- [ [ dqdp + 0{h l ~ n ) , 

{27Th) n J J\ p 2 +v{q) _ E \ <& 

see |14[ l2"H] and references given there. We note that this implies a precise upper bound 

(2.3) # { Spec(tf) n [E - Ch, E + C/i] } = ^(/i 1 -™) , 

which can be improved further by making assumptions on the classical flow of the Hamil- 
tonian p 2 + V(q) on Q, see [IHI25]. 

For open systems, with the simplest example given by the Hamiltonian in (|2.1)1 . real eigen- 
values are replaced by complex resonances. The simplest definition (easily made rigorous 
in the case (j2.1Jl ) comes from considering the meromorphic continuation of the resolvent. 
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Defining the Green's function G(z, q', q) for Imz > through 

{z-H)- l u{q')= ( G{z,q',q)u{q)dq, u E C c °°(M n ) , 

then G(z,q',q) admits a meromorphic continuation in z across the real axis. Its poles for 
Imz < (which do not depend on (q' , q)) are the quantum resonances of H. 

Counting of resonances is affected by the dynamical structure of the scatterer much more 
dramatically than counting of eigenvalues of closed systems. Since we are now counting 
points in the complex plane we need to make geometric choices dictated by dynamical and 
physical considerations. Here we consider scatterers and energies exhibiting a hyperbolic 
classical flow, and regions in the lower half-plane which lie at a distance proportional to 
h from the real axis. This choice is motivated as follows. Quantum mechanics interprets 
a resonance at z = E — 27 in terms of a metastable state, which decays proportionally to 
exp(—tj/h). Hence for 7 3> h the decay is so rapid that the state is invisible. On the 
other hand, for many chaotic scatterers there are no resonances with 7 <C h. One class for 
which this is known rigorously consists in the Laplacian on co-compact quotients HP/T, 
H = — /i 2 Ae"/r, when the dimension of the limit set satisfies S(T) < (n— l)/2. This follows 
from the work of Patterson and Sullivan — see the discussion below and [37j. 

After a complex deformation (see [53] and references given there) the long living quantum 
states should semiclassically concentrate on the set of phase space points which do not 
escape to infinity, that is on the trapped set Ke defined as follows: let H# be the Hamilton 
vector field of the Hamiltonian H(q,p) = p 2 /2 + V(q): 

n 

E h = ^2pjd qj - d qj V(q)d Pj . 
i=i 

Then 

Ke == T + (E) fl F-(E) , with the forward/backward trapped sets 

def — 

r±(E) = {pe^E : exptS^(p) y 4 00 > — ► °°} ■ 

Suppose that the flow generated by Eh is hyperbolic near Ke 1 for E' close to a non- 
degenerate energy E. That means that the field S# does not vanish on the energy surfaces 
E E , = {p 2 + V(q) = E'} C T*R n for E' w E, and that for p e E E > near K E >, 

T P Z E > = RE H {p) © E+(p) © E-{p) , 

S E / 3 p 1 — > E±(p) C TpT^E' is continuous, 
(2 ' 5) d(exptE H )(E±(p)) = E±(expt~ H (p)) , 

\\d(exptS H )(X)\\ < Ce ±xt \\X\\ , for all X G E±(p), Tt > 0. 
Weaker assumptions are possible — see [201 §5] and j^Hl §7]. 
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Typically, the set K E has a fractal structure and in the semiclassical estimates the 
Minkowski dimension naturally appears: 



We say that K E is of pure dimension if the supremum is attained. For simplicity of the 
presentation we assume that this is the case. 

Under these assumptions the estimate (|2.3|) has an analogue for chaotic open systems 
[HH] . For Co > there exists C\ such that 

(2.6) #{Res(H)f]{z : \z - E\ < C h } } < CJT** , dim K E = 2fi E + 1 • 

We notice that for a closed system the trapped set is the entire energy surface, so that 
in that case \i E = n — 1, hence ()2.6|) is consistent with ()2.3|) . In this note we use open 
quantum maps to provide the first evidence that this precise estimate is optimal. 

We should also mention that, as was already stressed in the work of Sjostrand j3U], the 
estimates involving the dimension are only reasonable when the flow is strictly hyperbolic. 
In the case of more complicated flows the estimates should be stated in terms of properties 
of escape or Lyapunov functions associated to the flow - see [HHllSSj- For expository reasons 
the estimates involving the dimension are however most persuasive. 

2.2. Survey of related results. The first indication that fractal dimensions enter into 
counting laws for quantum resonances of chaotic open systems appears in a result of 



where fh is any number greater than the dimension of the trapped set in the shell H~ 1 (E — 
I/C2, E + I/C2). In a homogeneous situation, such as for instance obstacle scattering, the 
dimension of K E , 2fi E + 1, is independent of E, so that m > 2(/j lE + 1). 

The improvement in quoted in ()2.6|) lies in providing a bound for the number of 
resonances in a smaller region D(E,Ch) = {z £ C : \z — E\ < Ch}. Heuristic arguments 
suggesting that the estimate ()2.7|) should be optimal were given in |3T] and [32J. 

Another class of Hamiltonians with chaotic classical flows and fractal trapped sets is 
given by Laplacians on convex co-compact quotients, M/T. Here T is a discrete subgroup 
of isometries of the hyperbolic plane H, such that 

• All elements 7 G T are hyperbolic, which means that their action on EI can be 
represented as 



dimK E = 2n — 1 — sup {c : limsupe c vol{p G T, E : dist(K E , p) < e} < 00} . 



Sjostrand [30] : 



(2.7) 



#{Res(#)n{,2 : \z-E\<5, Im z > -7 } } < C x 8 y 
Ch<j<l/C, max(h$,h/j) < 5 < 2/C , 




(2.8) 



ao 1 o a - 1 {x,y) = e e{l) {x,y), (x, y) G B. ~ R + x R , a G Aut(H) . 
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• If 7T : H — > H/r, and A(T) C is the limit set of T, that is the set of limit points 
of {7(2) : 7 G T}, 26I, then 7r(convex hull A(T)) is compact. 

The trapped set is determined by A(T): trapped trajectories are given by geodesies con- 
necting two points of A(r) at infinity, and 

dim K E = 25 T + 1 , S T = dim A(T) . 

The limit set is always of pure dimension, which coincides with its Hausdorff dimension. 
A nice feature of this model is the exact correspondence between the resonances of 

# = /i 2 (-A H/r -l/4), 

and the zeros of the Selberg zeta function, Z r (s) 2 : 

(2.9) zeRes(H) Z v {s) = , z = h 2 (s(l - s) - 1/4) , Res<S r , 

where the multiplicities of zeros and resonances agree. The Selberg zeta function is defined 
by the analytic continuation of 

{7} k>0 

where {7} denotes a conjugacy class of a primitive element 7 G T (an element which is 
not a power of another element), and we take a product over distinct primitive conjugacy 
classes (each of which corresponds to a primitive closed orbit). The length £(7) of the 
corresponding closed orbit appears in (J2.8|) . The exact analogue of (|2.6|) is given by 

(2.10) #{s:Z r (s) = 0, Res > -C , r < Im s < r + C x } < C 2 r Sr , 

which is a consequence of an estimate established by Guillope-Lin-Zworski j2D] i n a more 
general setting of convex co-compact Schottky groups in any dimension, 

(2.11) \Z T (s)\ < C K e CK ^ Sv , Res > -if, for any K. 

This improved earlier estimates of |59j . the proof of which was largely based on [oTI] . 

In the (non-quantum) context of rational maps on the complex plane, similar results were 
obtained concerning the zeros of associated zeta functions |11[ 134*] . Take / a uniformly 
expanding rational map on C (for instance z 1— > z 2 + c, c < —2), and call f n its n-fold 
composition. The zeta function associated with this map is given by 

(2.12) z(s) = ex P (-f>- £ jcrwr 



n=l f n (z)=z IW 



2 We refer to 02 f° r this and a general treatment. The term \ in the definition of the Hamiltonian H 
comes from requiring that the bottom of the spectrum of H is 0, so that Green's function (H — A 2 ) -1 is 
meromorphic in A G C 
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Then the number or resonances in a strip is also given by a law of the type (|2.10|) . where 
Sr is replaced by the dimension of the Julia set: 

J= \J{z : Hz) = z}. 

n>\ 

Note that this set is also made of "trapped orbits" . 

2.3. Survey of numerical results. The first model investigated numerically was perhaps 
the hardest to give definitive results. Lin [3U1 HH] studied the semiclassical Schrodinger 
Hamiltonian (|2.1|) with a potential made of 3 Gaussian "bumps". The semiclassical res- 
onances were computed using the method of complex scaling and were counted in boxes 
of type [E — 8, E + S] — i[0,h] with 5 fixed. The purpose was to verify optimality of 
Sjostrand's estimate ()2.7|) with these parameters. The results were encouraging but not 
conclusive. Since for small values of h the method of [3U] required the use of large matrices 
to discretize the Hamiltonians, the range of /i's was rather limited. 

A different point of view was taken by Lu-Sridhar-Zworski j22] where resonances for 
the three discs scatterer in the plane were computed using the semiclassical zeta function 
of Eckhardt-Cvitanovic, Gaspard, and others (see for instance ^21 CHI EE] and references 
therein). The zeta function is computed using the cycle expansion method loosely based 
on the Ruelle theory of dynamical zeta functions. Although it is not rigorously known if 
the resonances computed by this method approximate resonances of the Dirichlet Lapla- 
cian in the exterior of the discs, or even if the semiclassical zeta function has an analytic 
continuation, proceeding this way is widely accepted in the physics literature. Resonances 
z = h 2 k 2 were counted in regions 

(2.13) { k G C : 1 < Rek < r, Im£;>-7}, r -> oo , 

which under semiclassical rescaling correspond to counting in [1/2,2] — i[0, jh/2], h — ► 0. 
Let us denote the number of resonances (zeros of the semiclassical zeta function) in ()2.13|) 
by N(r,j). The fractal Weyl law corresponds to the claim that for 7 large enough, 

(2.14) N(r, 7) ~ C(rf) r^ +1 , r -> 00 , 

where 2/i + 1 is the dimension of the trapped set in the three dimensional energy shell (for 
such homogeneous problems, all energy shells are equivalent). In [32] the prediction ()2.14|) 
was tested by linear fitting of log N(r, 7) as a function of logr: 

log iV(r, 7) = (0(7) + l) logr + 0(1). 

We found that the coefficient 0(7) was independent of 7 for 7 large enough, and that it 
agreed with /1. The counting was done for three different equilateral disc configurations, 
parametrized by p = R/a where a is the radius of each disc, and R the distance between 
them. We also noticed that if 7 P is the classical rate of decay for the p configuration, then 

a p {x-f p /2) 

Hp 
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is essentially independent of p for 1 < x < 1.5. This corresponds to a numerical observation 
that for each p the distribution of resonance widths (imaginary parts) peaks near 7 = 7 P /2. 

Encouraged by the results of |3Z|, the cycle method was used in j20] to count the zeros of 
the Selberg zeta function for a certain Schottky quotient, but the results were not definitive. 
For the dynamical zeta function (J2.12j) with f(z) = z 2 + c, c < —2, the resonances were 
computed by Strain- Zworski |54J, using a different method based on the theory of the 
transfer operator on Hilbert spaces of holomorphic functions introduced in [20 . The zeros 
were counted in a region of the same type as in (|2.13p . 

{s : Res>-K, 0<lms<r}, 

where real parts and imaginary parts exchange their meaning due to different conventions 3 . 
By reaching very high values of r we saw a very good agreement of the log-log fit with the 
fractal Weyl, with /1 given by the dimension of the Julia set. 

In the model considered in this paper, we can verify the optimality of the fractal Weyl 
law on much smaller scales (see Table Q and the numerics presented in |4TH EHl ESI ) • That 
could not be seen in the other approaches. 

2.4. Related models in physics. The behaviour of quantum open systems has been re- 
cently investigated in situations where the classical dynamics has chaotic features. The 
physical motivation can originate from nuclear or atomic physics (study the lifetime sta- 
tistics of metastable states, possibly leading to ionization), mesoscopic physics (study the 
conductance, conductance fluctuations, shot noise in quantum dots or quantum wires), 
and from waveguides (optical wave propagation in an optical fiber with some dissipation, 
microwave propagation in an open microwave cavity). 

2.4.1. Kicked rotator with absorbing boundaries. In [31[7| a kicked rotator with absorbtion 
was used to model the process of ionization. The classical kicked rotator is Chirikov's 
standard map on the cylinder, which is a paradigmatic model for transitions from reg- 
ular to chaotic motion [§]. Quantizing the map on L 2 (T X ) results in a unitary opera- 
tor U, a first instance of quantum map. To model the ionization process which takes 
place at some threshold momentum p ion , the authors truncate the map U to the subspace 
^ion = span { \pj) : \pj\ < Pi on }: a particle reaching that threshold is ionized, or equiva- 
lently "escapes to infinity" . Here the discrete values pj = 2nhj are the eigenvalues of the 
momentum operator on L 2 (T X ); the space 7i ion is thus of dimension iV « p ion /nh. This 
projection leads to an open quantum map, namely the subunitary propagator U- lon = n ion ?7, 
where n ion is the orthogonal projector on Hi on . For the parameters used by [3], the classical 
dynamics is diffusive, meaning that a particle starting from p = will need many kicks to 
reach the ionization threshold. 



3 Although frustrating, the different conventions of semiclassical, obstacle, and hyperbolic scattering 
show how the same phenomenon appears in historically different fields. 
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The matrix {7 ion was numerically diagonalized for various values of h with p ion fixed, 
and the distribution of the N level widths 7$ = — 2 In | A« | , Aj G Spec(f/i on ) was found 
approximately independent of h, such that the number of resonances 

n(iy, 7 ) = #{ 7l < 7 } 

scales like C{^)N in this case. In subsequent works [HI EDI UZj > this distribution was shown 
to correspond to an ensemble of random subunitary random matrices, more precisely the 
ensemble formed by the [aN] x [atN] upper-left corner (0 < a < 1 fixed) of a large N x N 
matrix drawn in the Circular Unitary Ensemble (that is, the set U(N) equipped with Haar 
measure). 

2.4.2. Quasi-bound states in an open quantum map. Recently, Schomerus and Tworzydlo 
|4*U] have performed a similar study for the quantized kicked rotator on the torus (obtained 
from the map of the former section by periodizing the momentum variable). They also 
"opened" the map by assuming that particles reaching a certain position window g 6 I 
"escape to infinity". The quantum projector associated with these "escape window" is 
denoted by Ul, so that the remaining subunitary quantum map reads U op = (I — Hl)U. 
The main difference with the case studied in the previous section lies in the strongly chaotic 
motion (as opposed to diffusive), due to a different choice of parameters. The map has a 
positive Lyapunov exponent A, and a typical trajectory will escape after a few kicks: the 
average "dwell time", called td, is of order unity. 

The eigenmodes associated with eigenvalues bounded away from zero are called "quasi- 
bound states" , as opposed to the "instantaneous decay modes" associated with very small 
eigenvalues. The authors provide numerical and heuristic evidence that, in the semiclassical 
limit, the number of quasi-bound states grows like N cS = N 1 ~ 1 ^ Xtd \ This shows that 
most eigenvalues of U op are very close to zero, while only a small fraction N e g/N remains 
bounded away from zero. The authors also plot the distribution of the ~ N eS quasi-bound 
eigenvalues: again, it resembles the spectrum of a random subunitary matrix obtained by 
keeping the upper-right corner block of size iV e g- of a [T£>N e s] -dimensional random unitary 
matrix. 

The quantized baker's relation we will study in £ |4.6f J51 will be of similar nature. For the 
map (|4.38|) . the fractal dimension v given in (|3.7|) can be shown to be close to the formula 
1 — 1/(Atd), in the limit when the dwell time td is large compared to unity (limit of small 
opening). 

2.4.3. Conductance through an open chaotic cavity. The "scattering approach to semiclas- 
sical quantization" [3J H3J E], consists in quantizing the return map on a Poincare 
surface of section of the Hamiltonian system under study. Within this approach, the scat- 
tering matrix of the open system can be expressed as a "multiple-scattering expansion" in 
terms of the quantized return map. 



12 
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Using that framework, Beenakker et al. j^Zj study the quantum kicked rotator defined 
in the previous section, in order to understand the fluctuations of conductance through 
a quantum dot. The evolution inside the closed dot is represented by the same unitary 
matrix U as in last subsection, and its opening L is split into two intervals, L2 and L\, 
which represent the two "leads" bringing in and taking out the charge carriers from the 
dot. The orthogonal projector corresponding to these openings reads I1 L = U Ll ©Il i2 . The 
conductance can then be analyzed from the scattering matrix of the dot: 

(2.15) = n x {e-** - u{\ - n^}- 1 ^ . 

Here d G [0, 2ir) is called the quasi- energy. In terms of this parameter, the "physical half- 
plane" corresponds to Imi? > 0: the matrix 5($) has no singularity in this region. On the 
opposite, the resonances analyzed in the previous section, which are the poles of S("d), are 
situated in the region Imfi < 0. 

While is unitary, its subblock t = f H^S^Il^ describes the transmission from the 
lead L\ to the lead L2. The dimensionless conductance (which depends on is given by 
the Landauer-Biittiker formula g = tr (tt*). The eigenvalues of tt* (called "transmission 
eigenvalues") can be either close to 1 (corresponding to a total transmission), or close 
to (corresponding to a total reflection), or inbetween. The last case corresponds to 
genuinely quantum transmission eigenmodes, which are partly transmitted, partly reflected, 
due to interference phenomena inside the dot. The "quantum shot noise" is due to these 
intermediate transmission eigenvalues. A simple measure of that noise is given by the Fano 
factor [0] F = tr(££*(l — tt*)) / ti tt* . Using similar arguments as in the former section, the 
authors show that the number of intermediate transmission eigenvalues also scales like N e g, 
and thereby estimate the Fano factor, by assuming that these eigenvalues are distributed 
according to the prediction of random matrix theory. 

In £jf3we will analytically compute both the conductance and the Fano factor in the case 
of the open quantum relation B h . 



3. Classical dynamics 

3.1. Symplectic geometry on tori. We consider the simplest class of compact symplectic 
manifolds, the tori, 

n 

T 2n dcf E 2ry Z 2n ~ (fl x I) n , u = ^dq t f\ dp t , {q, p) G T 2n . 

t=l 

Here and in what follows, we identify the interval I = [0, 1) with the circle T 1 = R/Z. 
A Lagrangian (submanifold) A C T 2n is a n-dimensional embedded submanifold of T 2ra 
such that c<j|a = 0. We recall the following well known fact (see for instance Theorem 
21.3.2]): 
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Proposition 3.1. Suppose that A C T 2n is a Lagrangian submanifold, and that {qo,Po) G A. 
Then, after a possible permutation of indices, there exists k, < k < n, and a splitting of 
coordinates: 

q = (4, q") , v = (p',p") , q = (qi, ■■■qk), v" = (pk+i, ■■■ ,p n ), 

such that the map 

A3 (q,p) i— > (q",p') El n ~ k xl k 

is bijective from a neighbourhood V of (qo,po) to a neighbourhood W of (q'o,p'o). Conse- 
quently there exists a function, S = S(q",p') defined on W , such that A n V is generated 
by the function S , that is, 

AHV = { (d p ,S(q",p'),q";p',-d qll S(q",p')) , (q",p') G w) . 

In this paper we will also consider singular Lagrangian manifolds obtained by taking 
finite unions of Lagrangians with piecewise smooth boundaries. 

3.2. Symplectic relations. 

3.2.1. Symplectic maps. A symplectic (or "canonical") diffeomorphism on the torus T 2 " is 
a diffeomorphism k : T 2n — > T 2n which leaves invariant the symplectic form on T 2n : 
k*uj = to. An equivalent characterization of such a map is through its graph T, which is the 
2n-dimensional embedded submanifold of T 2n x T 2n , defined as 

T K = {(p';p) : p = (q,p)eT 2n , p' = «(p) } . 

Using the identification I" = M n /Z n , we setup the reflection map I n 3 p i— > —p G I™, and 
define the twisted graph O §25.2] 

(3.1) K = {{q',q;p',-p) : (q',p';q,p) G T K } C T 4n . 

Then the diffeomorphism k is symplectic iff T' K is a Lagrangian submanifold of T 4n (equipped 
with the symplectic form Y^l=i dq'j A dp'j + dqj A dpj). For this reason, we will sometimes 
denote T' K by A K . 

The definition of the twisted graph is clearly dependent on the choice of the splitting of 
variables {q,p), which will be related to a choice of polarization in the quantization process. 

More generally, one can consider invertible maps on T 2n which are smooth and symplectic 
except on a negligible set of singularities (say, discontinuities on a hypersurface) . The 
twisted graph of such a map is then a singular Lagrangian submanifold of T 4n . 

Example. The usual "baker's map" is the following piecewise-linear transformation k on 
T 2 : 



(3.2) K(<?,p) d = 



(2q,p/2) if < g < 1/2 

(2g-l,p/2 + l/2) ifl/2<g<l. 
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The twisted graph of k: 

A K = f { {q, q; p , -p) : (q,p) G T 2 , (q',p') = K(q,p) } 

is a singular Lagrangian submanifold of T 4 . It can be decomposed into A K — Ao U Ai, with 
the components 

A i =|(2g-j,g;^ > -p) : j/2 < q < j/2 + 1/2, p G 1 1 

= {(2g-j,g;p',-2p , +j) : j/2 < < j/2 + 1/2 } . 
ngian in T 4 

a 2-dimensional square. 



Each Aj is locally Lagrangian in T 4 and, as a manifold with corners, it is diffeomorphic to 



3.2.2. Multivalued symplectic maps. A canonical (or symplectic) relation is an arbitrary 
subset T C T 2n x T 2n , such that 

r' = {(q',q; p',-p) ■ (q',p';q,p) er} 

is a Lagrangian submanifold of T 4n . 

We are interested in symplectic relations coming from multivalued symplectic maps. A 
multivalued map is the union of finitely many components Kj, where canonical 
diffeomorphism Kj between an open subset Sj with piecewise smooth boundary of T 2n and 
its image <S'- = Kj(Sj) G T 2n . A priori, the sets Sj (respectively <S'-) can overlap, and their 
union can be a proper subset of T 2n . 

Each map Kj is associated to its graph 

Tj = {(Kj(p);p) : peS 3 } , 

and the symplectic relation can now be defined through its graph 

r = U r ;' 

3 

or equivalently its twisted graph V (defined from Y as in (|3.1Jl ). V is a singular Lagrangian 
in T 4n . 

The inverse relation can be defined by 

r- ld ={(p;p') : (p';p)eT} = \J{(Kj 1 (p);p) : p G S'j } , 

i 

and the composition of two relations by 

f or d ^ f {(/;p)GT 4 " : 3 p' G T 2n , ^) G r and (p";p')ef}. 

Following j^J Theorem 21.2.4], we note that Tor will be a (locally) smooth symplectic 
relation if f x T C T An x T 4n intersects 

W,f/,pf,p) : p",p', P GT 2 ncT 4 ™xT 4 " 
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cleanly, that is the intersections of tangent spaces are the tangent spaces of intersections. 

We can then iterate a relation T, defining a multivalued dynamical system {T n : n e Z} 
on T 2n . In §3.41 we will give a stochastic interpretation to this system. 

3.3. Open baker's relation. The dynamics we will consider takes place on the 2-torus 
phase space, 

t 2 = {p = (i,p) 9, p e 1} . 

On this phase space, we define two vertical strips Sj (j = 1,2) from the data of four real 
numbers Di, D 2 > 1 and £ 1; l 2 > 0: 

(3.3) Sj = {{q,p) : qelj, pel}, with Ij = (-^, j = 1,2. 

The strips are assumed to be disjoint, which is the case if we impose the conditions: 

4 + 1^4 , h + 1 . 

< — - and — - — < 1 . 



D x ~ D 2 D 



2 



The corresponding baker's relation is made of two components Bj, j = 1, 2 associated with 
linear symplectic maps defined on the two strips: 

(3.4) Bj = { (//; p) : {q',p') = fa - Ij, P -^j ,p=(q,p)e Sj } . 

The baker's relation is defined as the graph B = B\ U B 2 . One clearly notices that each 
component map is a hyperbolic diffeomorphism, with positive stretching exponent log-Di 
(resp. logD 2 ). At all points where the map is defined, the unstable (stable) direction is 
the horizontal (vertical) one. 

Since the two strips are disjoint, each point p e T 2 has at most one image. In the 
notations of Proposition 13.11 (taking q" — q, q' — q'), each Lagrangian component Bj can 
be generated by the function 

(3.5) Sj(q,p') = Dj ^q — -^-J ^p' — -^-J defined on the square { (q,p') e Ij x Ij } . 
Let 

ir L , n R : fxf — > T 2 

be the projections on the left and right factors respectively. From the definition (|3.4|) . the 
set hr{B) = Si U S 2 is made of points on p e T 2 which have an image through the relation 
B. Hence, a point p e" ttr(B) is said to escape from the torus at time 1. Similarly, a point 
p e~ n L(B) = tir^B^ 1 ) is said to escape from T 2 at time —1. This "escape" is the reason 
why we call this relation an "open" relation: the system is not "closed" because it sends 
particles "to infinity", both in the future and in the past. 
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We define 

oo 

(3.6) T ± d ^ f fl tt r (BF») 

71=1 

the set of points which do never escape from T 2 in the past, respectively in the future. One 
checks that these subsets have the form 

r_ = C x I, r+ = I x C, 

where C C I is a "cookie-cutter set" in the sense of ^Bj : if we consider the two contracting 
maps on I 



Q 

IjW = 

this closed set is defined as 



f j (q) = ±— qel, J = 1,2 

U 3 



C = {q E I : f h o ■ ■ ■ o f jn (q) = q for some sequence j m E {1,2}}. 

The Hausdorff dimension of C (which is equal to its Minkowski and box-counting dimen- 
sions) is given by the unique < v < 1 solving 

(3.7) D\ v + D 2 U = 1 . 

The trapped set (or set of nonwandering points) is defined as the set of points which never 
escape from T: 

K = T + nT_=C xC, dimK = 2z/. 

The baker's relation is a hyperbolic invertible map on the set K , which is a "fractal repeller" . 
This relation is a model of Smale's horseshoe mechanism. 

The simplest case consists in considering a symmetric baker's relation, with D\ = D2 = 
D,£ = £ 1 = D-£ 2 -l- 

^ 3 ' 8 ^ i £+1 ( i> — ( — I 

^<l-?<^ =► (q', P ')=[D( q -l)+e+l, P ± +1 

Now C C I is a symmetric 1/D— Cantor set. Notice that if we take D = 2, £\ = 0, 
£ 2 = 1, we obtain the usual (closed) baker's map described in the example of Section I3~T1 
for which the trapped set (= T 2 ) has dimension 2. The "3-baker" relation described in 
(|1.3|) corresponds to D = 3, £ = 0. 

For such a symmetric baker's relation, the analog of the fractal exponent of (J2.6)) is: 

log 2 

He < ► v = 1 p: • 

logD 
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3.4. Weighted symplectic relations. To give a multivalued map T a physical meaning, 
we assign Markovian weights Pj{p) to the different "jumps", p h- > Kj(p). The associated 
dynamical system is then stochastic, each point p having finitely many images with well- 
prescribed transition probabilities Pj{p)- The sum of all the probabilities from p must 

satisfy < P{p) == J2j Pj(p) — 1> so that (1 — P(p)) is the probability that p "escapes to 
infinity" . 

The weights associated with the inverse relation T _1 are the same: each point p' jumps 
back to kJ (/?') with probability -Pj(p') = Pj(n~ 1 (p')). Hence, the weights must also satisfy 

Such a weighted relation (in geometric optics one would speak of a "ray-splitting" map) 
induces a discrete-time evolution of "mass distributions", which is in general dissipative: 
the full mass decrease at each step, the system expelling part of the mass "to infinity" . 

In more mathematical terms, we assume that the symplectic relation T C T 2n x T 2n 
comes with a nonnegative measure (or weight) p on T, which for any Xa G C°°(T 2n , [0, 1]), 
a = L,R, satisfies 

(3.9) ^lXL^ R XRP) = gl LXR -r, g* LXR eC™(T 2n ), 0<g?**<l, 

n! 

where ttl, 7Tr : V — > T 2n are projections on left and right factors respectively, and uo is 
the symplectic form on T 2n . The condition (J3.9)) implies that 7r Q |r is a local bijection, 
which forces T to be a piecewise smooth union of graphs of symplectic transformations, as 
defined in ^3.2.21 When V is singular, that is a union of smooth symplectic relations with 
boundaries, we demand that 

g xLxn eC o 0( j2^ if SU pp(7r2xr^ X Ji)nar = 0, 

where dT is the union of the boundaries of the smooth components. 

The reason for introducing the measure p is to have a quantity independent of the choice 
of coordinates on T. On T 2n , an obvious intrinsic measure is given by the symplectic form, 
hence g^ LXR are well defined. Building an atlas of the manifold V we can use these functions 
to describe p in local coordinates. 

We denote a weighted relation by (T, p) As explained above, one can invert such a 
relation, as well as compose them. 

If (/?'; p) G r \ dT, the probability of a transition from p to p' — K^{p) is obtained by 
letting xr (resp. xl) be supported in a sufficiently small neighbourhood of p (resp. of p'), 
with Xr(p) — 1> Xl(p') = 1- This probability is then given by 

(3.10) P n (p) = gx>**(p) = gl^(p') = P' h {p') . 

Examples. The simplest example is given by a graph of a symplectic transformation k : 
T 2n — > T 2n in which case the density p is obtained by taking p = n* L (uj n /n^ = ir^ t (uj n /n\), 



18 S. NONNENMACHER AND M. ZWORSKI 

where the equality follows from k*u = u. A slightly more complicated example is given by 
taking a union of two non- intersecting graphs Tj of Kj, j = 1, 2, and putting 

H = (7c R \ ri y(g 1 u n /n\) + (7c R \ T2 y(g 2 u n /n\) , 

where gj G C°°(T 2n ; [0, 1]) satisfy g% + g 2 < 1 and g\ o k^ 1 + g 2 o k 2 1 < 1. In this case, 
9j(p) = Pj(p)- 

In the case of an open baker B defined in £ J3.3[ for instance the symmetric 3-baker (|1.3j) . 
a natural p comes from pulling back the Liouville measure uo to each component Bj given 
in (|3.4j) . One obtains 

(3.11) TTii* /i = l/iu/ 2 (g) dp , n L *p = l huh (p') dq dp' . 

These equations fully determine the measure p on B. 

A more interesting example, which will be relevant in is given by the following 
multivalued generalization of the symmetric 3-baker: 

2 2 2 

B = U(5 + (0^/3;0 I 0)) = [JU%. where 

(3.12) 1=0 k=ij=o 

Bkj = |(3g,^;g,p) : q e 4, p g I J , h = (0,1/3), J 2 = (2/3,1). 

Each point p G <Si U 1S2 = (I\ U I 2 ) x I has 3 images, and each point p' G T 2 has two 
preimages. 

The following measure on B will arise in the quantum model studied in We define it 
explicitely on each component B^, using the right projection on S^'- 

-1 sm 2 (vrp) 

J 9sm (vr(p + j)/3) 
sin (7rp) 

n R*HB 2l = o ■ 2/ / 7 ■ ^n^j ? rf 9 OP, J = 0,1,2. 

2j 9 sin (7r(p + j — 2)/3) 

The functions on the right hand sides are the probabilities Pj(p). The sum of these com- 
ponents reads 

Kr * » = ( H ■ 2 T/JT -/ox ) ^lUIa (?) rf <? = 1/iUJa (?) d<? rf P • 

\ 9 7^ sm ^(p/s + J/ 3 )/ 

Here we used the fact 4 that sin 2 (_Dx)/ sin 2 (a; + jn/D) = D 2 , with D = 3 and 

x = np/3. This right pushforward is identical to that of ([3.11)1 : in both cases, any point 
p G (S\ U S 2 ) has an empty escape probability, 1 — P(p) — 0. 

4 The value of the sum at x = is equal to D 2 , and the sum is invariant under translation x 1— > x + kir/D. 
Fejer's formula for the Cesaro mean of the Fourier series shows that the sum is a trigonometric polynomial 
of degree D — 1 in x, hence it is constant. 



DISTRIBUTION OF RESONANCES FOR OPEN QUANTUM MAPS 19 

On the opposite, the left pushforward of p is given by 

sin 2 3717/ f 1 1 \ . . 

ttl* fi = — — - + • 2 / , — 7777 rf <? ^ • 

9 \sm rep' sm 7r(p' — 2/3) / 

Almost any point p' G T 2 has a nonzero escape probability through This left push- 

forward is obviously different from that of p. 

4. Quantized maps and relations 

Before giving the definition of the quantized baker's relation, we need to define the 
quantum Hilbert space corresponding to T 2 , as well as the algebra of quantum observables. 

4.1. Quantized tori. The quantization of tori T 2n = M 2n /Z 2 " has a long tradition in 
mathematical physics j2H H31 EJ- It can be considered as a special case of the Berezin- 
Toeplitz quantization of compact symplectic Kahler manifolds — see j22j and references 
given there. Here we will give a self-contained presentation of the simplest case from the 
point of view of pseudodifferential operators. 

We first recall from [T3j the quantization of functions / G C^(T*R n ), 

C£°(T*R n ) = {/ e C°°(T*R n ) : Va, /3 G N", sup \d^d^f{q,p)\ < 00} . 

(g,p)eT*R n 

To any / G S(T*~R n ) we associate its h-Weyl quantization, that is the operator f w (q, hD) 
acting as follows on %p G 5(M n ): 

(4.1) [r ^ hD )^ {q) ^ -L- J J f^+L^yi^r, P )^ r)drdp 

This operator clearly has the mapping properties 

f w {q,hD) : S{R n ) — ► S{R n ) , f w {q,hD) : S'{R n ) — > S'{R n ) . 

It can be shown [T4"| Lemma 7.8] that / 1— > f w (q, hD) can be extended to any / G C^(T*R n ), 
and that the resulting operator has the same mapping properties. Furthermore, f w (q, hD) 
is a bounded operator on L 2 (R n ). 

We now introduce quantum spaces associated with the torus T 2n . For that aim, we fix 
our notations for the semiclassical Fourier transform on iS'(IR n ): 

and as usual in quantum mechanics, Th^ip) is the "momentum representation" of the state 
ip. The torus quantum space is made of distributions ip G «S'(R n ) which are both periodic 
in position and momentum: 

(4.2) + =V(g), F h tP(p + £) = F h 4j(p). 

Let us denote by Ti 1 ^ this space of distributions. We have the following elementary 
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Lemma 4.1. Ti 1 ^ ^ {0} if and only if h — (2-kN)^ 1 for some positive integer N, in which 
case dim7i^ = N n and TC^ is generated by the following basis: 



(4.3) H n h = span { -= V%-f- j/N) : j G (Z/NZ) n 



The distributions elements of this basis will be denoted by 

j 

(4.4) \Qj) , Qj — — G I n is the position on which that state is microlocalized. 

One can check that for such a value of h, the Fourier transform Th, maps H. 7 ^ to itself. In 
the above basis, it is represented by the discrete Fourier transform 

— 2m(j j') IN 

(4-5) {FN)ir = e Nn/2 , J,fe(Z/NZr. 

It is also easy to check the following 

Lemma 4.2. Suppose that f G C£° (IR n x IR n ) satisfies f(q + £,p + m) = f(q,p) for any 
£, m G Z n . TTien t/ie operator f w (q, hD) maps ?{% to itself. 

Identifying a function / G C°°(T 2n ) with a periodic function on IR 2n , we will write Op h (f) 
for the restriction of f w (q, hD) on 7i^, 

c oo (T 2n )9/ Op fc (/) G £(WJf) . 

We remark that Op /l (l) = Id. The vector space can be equipped with a natural Hilbert 
structure. 

Lemma 4.3. There exists a unique (up to a multiplicative constant) Hilbert structure on 
H n h for which all Op h (f) : H n h H n h with f G C°°(T 2n ; R) are self-adjoint. 

One can choose the constant such that the basis in (14.3(1 is orthonormal. This implies 
that the Fourier transform on Ti 1 ^ (represented by the unitary matrix (14.5(1 ) is unitary. 

Proof. Let (•, •)o be the inner product for which the basis in ()4.3(1 is orthonormal. We 
write the operator f w (q, hD) on ?{% explicitely in that basis using the Fourier expansion of 
its symbol: 

For that let Le tTn (q,p) = (£,q) + (m,p), so that 

f w (q,hD)= ]T f(e,m)ex V (2niLl m (q,hD)). 

Applying this operator to the distributions ((4.4(1 . we get 

exp (2mLl m (q,hD)) \Qj) = exp {^(2{j,£) - (m,£))) |Q,-_ m ) , 



DISTRIBUTION OF RESONANCES FOR OPEN QUANTUM MAPS 



21 



and consequently, 

r(q,hD)\Q j )= F mj\Q m ), 

m&Z n /(NZ) n 

Fmj= E RtJ-rn-rN)(-l)M exp {j + m, £)) . 

Since 

F jm = f(-e,j-m + rN)(-l)^ exp ( - ^ (j + m, £>) 

= Yl fei-"»-^)(-l) M exp(^0- + m,l)), 

we see that for real /,/ = /, P/ m — -^mj- This means that f w (q,hD) is self-adjoint for 
the inner product (•, •) . We also see that the map / i— ► (-P 1 jm)i,m6(z/ivz) n is onto, from 
C°°(T 2n ; R) to the space of Hermitian matrices. 

Any other metric on TC% could be written as (u,v) = (Bu,v) = (u,Bv) . If (f w u,v) = 
(u, f w v) for all /'s, then Bf w = f w B for all /'s, and hence for all Hermitian matrices. 
That shows that B = eld, as claimed. □ 

This choice of normalization (•, •) can be obtained in a natural way, if we use the 
following periodization operator to construct ?{% from iS(M n ) [3]: 

Lemma 4.4. For any h = (27TN)" 1 , the periodization operator Pj 2 ™ : <S(R n ) — > defined 
below is surjective: 

(4.6) W>£<S(R"), [p T2n ^](g.) d J: f -L-^^--^), je{Z/NZ) n . 

In the rest of this article we will always assume that h = (2ttN)~ 1 for some N £ N, so 
the semiclassical limit corresponds to N —>■ oo. The scalar product on Ti 1 ^ will be (•, •) . 
From now on we will omit the subscript 0, and also often use Dirac's notation (»|«) for 
this product. For instance, the j-th component of a state ip £ in the basis f)4.4j) will 
be denoted by ip{Qj) = (Qj\ip). The Hilbert norm associated with (•, •) will simply be 
written || • ||. 

4.2. Lagrangian states. We want to characterize the semiclassical localization in phase 
space of sequences of states of the form ip = {ip h £ TCfyh^o- In general we will assume that 
each element of this sequence is normalized, \\iph\\ — 1> but all definitions can be extended 
to sequences such that the norms satisfy \\iph\\ — 0(h K ) as h — > 0, for some fixed K £ R 
(the sequence ip is then said to be tempered). 
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The localization of this sequence is first characterized through its microsupport, or wave 
front set, which is the following subset of T 2n : 

(4.7) WF fc (V>) = C { p G T 2 " : 3 / G C°°(T 2n ) , f(p) ? , || Op^J = 0(h°°) } , 

where C stands for the set theoretical complement. It is not hard to show Proposi- 
tion IV-8'] that this definition is equivalent to the following: p G" WF ^(ip) if and only if 
there exists a neighbourhood W p of p such that, for any / G C°°(T 2n ) supported in W p , 
II Op^(/)'0/i|| = 0(h°°). This yields the following 

Lemma 4.5. For any function f G C°°(T 2n ) with f = in an open neighbourhood of 
WF h (ip), we have || Op h (f)iph\\ = 0{h°°). As a consequence, the microsupport of a sequence 
ip = {4>h}, \\4>h\\ x h K , cannot be empty. 

Proof. The (compact) support of / can be covered by finitely many W Pi , and using a 
partition of unity associated with these sets we can decompose it as / = with 
supp(/j) C W Pi . We get the result by linearity, and using the second definition of WF^ (?/>). 

□ 

We also make the following observation: 

Lemma 4.6. Let ip = {i/j^ G 7i^}^o be a tempered sequence. Considering iph as a N n - 
component vector in the basis ([4.3)1 . we define ip h as the vector with complex conjugate 
components. Then 

WF h (j>) = {(q,-p) : (q, p) G WF h (ip)} . 
Proof. The definition (j4.1|) of Weyl's quantization gives, for any function / G C°°(T 2n ), 

Op h (f) $ = f W (Q, hD) 4, = f"(q,-hD)il> . 
The lemma follows from the definition (|4.7J) of the wave front set. □ 

Now let A C T 2n be a union of Lagrangian submanifolds of T 2n with piecewise smooth 
boundaries. 

Definition 4.7. A sequence of states ip = {iph G H^} is a Lagrangian state associated to 
A, which we denote by ip G 1(A), if for any M G N and any sequence of functions, 

G C°°(T 2n ) , l<j<M, /,|a = 0, 

we have 

(4.8) || o Ph (f M ) o.-.o opfcC/o m = o(h M ) Hvy . 

From the definition (|4.7jl of the microsupport, we obtain that, if the sequence ip is 
tempered, then 

(4.9) ipel(A) WF^)cA. 
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Indeed, suppose that p A. Then there exists / G C°°(T 2n ) such that / |a = and / = 1 in 
a neighbourhood of p. We can also find a G C°°(T 2n ) such that / = 1 on a neighbourhood 
of the support of a, and a(p) ^ 0. The symbol calculus (see [HI Chapter 7]) shows that 
for any M, Op h {a) Op h (f) M = Op h {a) + M (h°°). On the other hand || 0^ h (f) M ip h \\ = 
C(/i M ||?/'/ l ||), and as M is arbitrary and if) tempered, || Op h (a) ip h \\ = 0(h°°). In view of 
(gZZD, this gives (JOJ). 

We stress that the opposite implication in (|4.9J1 is not true in general. To see that 
consider n = 1 and the Lagrangian A = {(0,p) '■ p G 1} C T 2 . Let iph G 7i\ be the "torus 
coherent state at the origin" : 

MQi) = (^) 1/4 ^expi-nNiQj - r) 2 } , j = 0, . . . , N - 1 . 

Then one can check that ||^|| -^-A 1, that WF h (ip) = { (0, 0) } C A. On the other hand, 

|| Op h (sin (2Trq))if) h \ \ ~ ix\p2h , 

which shows that iph ^ -^(^-)- 

In the physics literature, Lagrangian states are usually called WKB states, and are intro- 
duced as Ansatze for eigenstates of integrable systems, using Bohr-Sommerfeld quantization 
formulae [2E]. For instance, in the case n = 1, if A is generated by the function S G C°°(T): 

(4.10) A s = {(q,-S'(q)), gel}, 
then for any function a(q) G C°°(I), the state iph £ defined as 

(4.11) AiQj) = ^-r=- exp(-2inNS(Qj)) , j = 0, . . . , JV - 1 , 

V JM 

is in I (As). In the next proposition, we generalize this construction to any dimension. 

Proposition 4.8. Let A C T 2n be an embedded Lagrangian manifold. Then for any p G A 
there exist Lagrangian states ip G 1(A), such that po £ WF^ (?/>). 

Proof. We take p = (9o?Po) £ A, and assume that there exists a neighbourhood V of p > 
and a function S G C°°(7r(V)) (where n(q,p) = q), such that A fl V = {(q; —d q S(q)) , g G 
7r(V)}. This is a particular case of Proposition 13.11 The general case of a generating 
function S(q",p') can be transformed to that of S = S(q) using the symplectic rotation 
(q',p ! ) i — ► (~p',q')- On the quantum mechanical side, this rotation is performed through 
a partial Fourier transform in the variable q'. Our construction below can be transposed 
to this general case through this Fourier transform (which acts covariantly on the Weyl 
quantization). 
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We also assume that the neighbourhood V is contained in the interior of I 2n , and we 
identify 7r(V) with a subset of I n . We first construct a Lagrangian state in L 2 (R™): 

(4.12) u h (q) = a(q)e-i s ®, 

with a symbol a G C°°(IR n ) compactly supported inside 7r(V), and such that a(q ) ^ 0. 
This state admits the norm ||w/i||z,2 = ||a||£,a. For any / G C°°(T 2n ), we apply the operator 
f w (q, hD) to that state. Although we could do it directly using (j4.1j) . we prefer to reduce 
the problem to the case of S = by conjugation with the unitary multiplication operator 

(4.13) v(q) i— > [e^ sw ^v\{q] = e^v(q) , 
where we can assume that S G C£°(IR n ). We then apply the operator 

G w (q, hD) d ^ ei^to) r(q, hD) e"*^) , 

to the function a(q). The symbol calculus shows that G(q,p) admits an /i-expansion, 
with principal symbol g(q,p) = f(q,P + d g S(q)): if / vanishes on A, then g vanishes on 
{{q,0) : q G ir(V)}. We get 

(q, hD) u h }(q) = G w (q, hD) a(q) = e'^W g"( q , hD) a(q) + 0(h) , 

The explicit integral 

[ g ^(q,hD)a}(q) = -^- J J g(^ L ,p) o(r) e*<™> drdp 

can be evaluated through the stationary phase method. The derivative of the phase vanishes 
at r = q, p = 0, so the integral admits the following expansion [231 §7.7] for q G 7i(V): 

(4.14) [g w (q, hD) a](q) = L (g a)(q) + hL x {g a)(q) + 0(h 2 ) . 

Here each function Lj(ga) is obtained by applying a certain differential operator (in (r,p)) 
on the function g((q + r)/2,p) a(r), taking the output at the point (r = q, p = 0). The 
first term is simply L (ga)(q) = g(q, 0) a(q). For q outside 7r(V), the nonstationary phase 
estimates show that 

(4.15) n ,, ftD ) M ,) = ((_A_)"). 

If f(po) 7^ 0, then Lo(g a) is nonzero in a neighbourhood W of qo, and we obtain 

(4.16) \\r(q,hD)u h \\ L 2 iMn) = \\g w (q,hD)a\\ LHRn) + 0(h) > \\L (g a) \\ L 2 {W) + O(h) . 

The left hand side is thus bounded from below by a positive constant. 

On the opposite, if / vanishes on A, then at each point q G 7r(V) we get L (ga)(q) = 0, 
which implies that \\f w (q,hD)uh\\L 2 (R n ) — 0(h). The same procedure can be iterated to 
show that, for any family of functions /j G C°°(T 2n ) vanishing on A, the function 

(4-17) u[ M \q) d = f h- M [fl o ■ ■ ■ o f?u h \{q) , 
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is uniformly bounded and smooth on M. n , and very small outside 7r(V), as in (|4.15|) . As a 
result, 

(4.18) m(q, hD) o • • • o /*(g, ^K|U W = /* M ||4 M) || = 0(/* M ) . 

We can now carry over the estimates (|4.16J4.18J) onto the state iph = Pj2n Uh G 7i 
where Px 2 ™ is the periodizing operator (|4.fij) . Since a(g) was supported inside ir(V) C I 
this state admits the following representation, which generalizes (|4.1U|) : 

(4.19) MQj) = = eM-^NSiQj)) , j G (Z/JVZ)" . 
The norm of this state is therefore the sum 



n 
h ■ 



HQ 3 )?= j dq\a(q)\ 2 + 0(h c 

[Z/NZ) n 



\m\ 2 = 

where we used the smoothness of a(q). Similarly, the projection on of the function 
defined in (I4.17J) satisfies 



(4.20) p^ u M( Qj) = ?h_^ + 0{l ^ h je(Z , 

From the smoothness of , we obtain the "projected version" of ()4.18|) : 

II Op h (f M ) o • ■ ■ o op^/o^ll = h M \\p^ n 4 A/) || = /i M ||4 M) IUw + o(h°°) = o(h M ) . 

On the other hand, if f(po) ^ 0, one easily deduces from (|4.16|) that 

l|Op fc (/)^|| = \\r{q,hD)u h \\ L 2 {un) + 0{h 00 ) >C + 0(h), C>0. 

These estimates show that the family ip G 1(A), and that p G WF^). □ 

Remark 4.1. The definition of 1(A) mimicks the Hormander-Melrose definition of La- 
grangian distributions [2U Definition 25.1.1] (see pQ for an adaptation to the standard 
semiclassical setting). The requirement that A is Lagrangian reflects the uncertainty prin- 
ciple, in the following sense. A Lagrangian submanifold is the lowest dimensional subman- 
ifold for which the conclusion of Proposition 14.81 holds, that is, for any p G A, there exists 
a state if) satisfying tp G 1(A) and p G WF/ t ( , 0). 

Indeed, let A be an embedded submanifold of T 2n . Let us assume that ip G 1(A), so (j4.8|) 
must hold for any family of functions fj\\ = 0. From the identity 

£[Op h (/,), OpfcC/,-)] = Op h ({fi, fj}) + 0(h) . 

we see that || Op h ({fi, fj}) ip h \\ = 0(h). As in the proof of ()4.9|) . we can show that if 
{fi, fj}(p) 7^ for some p G A, then p G" WFh(ip). Hence, if we want the conclusion of 
Proposition 14.81 to hold for A, then this submanifold must satisfy 

Wi, fj G C°°(T 2n ), /.U, £| a = {/i,/ i }| A = 0. 
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This property means that A is co-isotropic, and must be of dimension > n. Lagrangian 
manifolds are co-isotropic manifolds of minimal dimension. 

4.2.1. Singular Lagrangian states. We now give an example where A is a union of La- 
grangians with piecewise smooth boundaries (in fcj.2l we called such A a singular La- 
grangian). Let A s be given by (|4.1(Jj) and iph by (j4.11|) . Let us truncate iph to some 
proper subinterval [Q, Q'} C I, that is, replace the symbol a(q) by the discontinuous func- 
tion d(q) = a(q)liQ^Q>] (q). That gives a state iph E Ti-\- One could expect iph to be a 
Lagrangian state in I (As) (as is iph), or rather in /(As), where 

A 5 d ^ f A s n([Q,g , ]xi). 

This is not the case: one needs to include in the Lagrangian the singularity set 

A shlg = {(Q, P ) : P eI}u{(Q',p) -.pel}, 

which is the "periodized" conormal bundle of the boundary dA$. We will indeed prove 
that iph E I(A S U A sing ), which can be considered as a semiclassical, discrete analogue of 
singular Lagrangian distributions of Guillemin-Uhlmann [HJ] and Melrose-Uhlmann |34| . 
We have the following 

Lemma 4.9. Let us truncate the state (I4.19J) to a hypercube H C I n , H = YYi =1 [ote, 

(4.21) iPhiQj) = am N l% m expi-XnNSiQj)) , j E (Z/iVZ)" . 

Then iph is associated with the singular Lagrangian A s UA sing , where A s = {(q, —d q S(q)) , q E 
H} and 

n 

(4.22) A sing = (J ({(q,p) : q e = a e , p t El, q m E [a m ,(3 m ], p m = -d qm S(q), m ^ t) 

U {(q,p) : qi = fa, Pe El, q m E {a m ,(3 m }, p m = -d qm S(q), m ^ £}) . 

Remark. It would be tempting to generalize the lemma by replacing the hypercube H by 
an arbitrary set S with smooth boundaries. However, if n — 2, S = 0, a = 1, and dS does 
not contain a segment with rational slopes then 

WF h {i> h ) = [S x {0}) U (dS x I 2 ) . 

The second component being 3-dimensional, this set is certainly not contained in a finite 
union of Lagrangians. 

Proof. As in the proof of Proposition I4.8[ we can, by conjugation with the operator ()4.13|) . 
reduce the proof to the case S = 0. We first consider states defined on W 1 , localized on the 
hypercube H C M n : 

(4.23) u h (q) = l H (q)a(q), a E C 00 (R™) . 
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We use the following 

Lemma 4.10. Let A = H x {0} and A sing be as in Lemma \4-9\ The ideal J of periodic 
functions vanishing on the singular Lagrangian A U A sing is (infinitely) generated by 

9jiP, q) = s in {ir(qj - sin (irfa - /%)) sin(7rp J ) , 

9ij{Q,p) = sin(Trpi) Bm(irpj) i^j, 1 < i, J < n , 

<t>M,v) = <P(.Qj,Pi, ■ ■ ■ ,Pj-i,Pj+i, ■•■ ,Pn), where <f>(,Qj, •) = 0, OLj < qj < (3j, 
ip(q) , where ip G C°°(I n ) vanishes on H . 

Proof. We only give the proof for the following model (n — 2), which contains all the basic 
ingredients of the general case. Let us study the ideal of functions vanishing on 

(4.24) =p 2 = 0}U {q 2 = pa = 0} U { Pl = p 2 = 0}) n {q x > , q 2 > 0} . 

The functions vanishing on the first factor in the intersection are generated by qipi , q 2 p 2 , 
and p\p 2 . Writing an arbitrary function F(q,p) as 

F(qi,q 2 ,Pi,p 2 ) = F (p u p 2 ) + qiF^, q 2 ,p 2 ) + q 2 F 2 (q 1 ,q 2 ,pi) + 

+ qiPiF n (qi,q2,Pi,P2) + q2P2F 22 (q 1 ,q 2 ,p ll p 2 ) , 

we need to find conditions for qiF 1 (qi,q 2 ,p 2 ) and q 2 F 2 {qi, q 2 ,Pi) to vanish on ([4.24)1 . We 
treat the first function by expanding it as 

Fi{qi,q 2 ,P2) = F 10 (q 1 ,q 2 ) + p 2 F 12 (q 1 ,p 2 ) +^2-^122(91,92,^2)- 

This forces F w (qi,q 2 ) to vanish identically in {qi,q 2 > 0} and Fi 2 (qi,p 2 ) to vanish identi- 
cally in {qi > 0}. 

The function F 2 (qi, q 2 , Pl ) is treated identically. Hence the functions vanishing on (|4.24j) 
are generated by qiPi, q2P2, P1P2, and all the smooth fuctions ip(qi, q 2 ), (pi(qi,p 2 ), $2^2, Pi) 
vanishing on {qi, q 2 > 0}. The transposition to the torus setting gives the lemma for that 
case. The general case can be proven similarly. □ 

This lemma means that any F G J can be decomposed as 

F = J2 hi 9a + 9j + fj Oj + V) , 

where the functions /. are smooth and either periodic or antiperiodic in each variable, so 
that f,g, are periodic in all variables. 

The action of each term (/ g) w (q, hD) on the state (J4.23|) can be written 

(fg) w u h =((f a) w o g w + h L(f, a, g)) t H , 

where L(f,a,g) is a pseudodifferential operator of norm 0(1). Therefore, we are reduced 
to study the action of the generators g w (q, hD), g = g^ , gj , 4>j,ip, on the characteristic 
function 1h(q)- 
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We first note that ip 1# = (pjln = 0, so there is nothing to prove in this case. 

For each j e { 1, . . . , n }, the generator g<j contains a factor sin(7rpj). Up to an error 
0(h), we first quantize this factor and apply it to H#: 

1 dof 

sin(irhDj)l H {q) = — (l H (qj + nh,q') - l H (qj - Kh,q')) = bj(q) . 

The function bj(q) is supported in the strips Sj = { \ qj — <x,| < irh } U { \qj — /3j\ < irh }, 
where it takes values ±1. We now apply the remaining factors of gf this amounts to 
multiplying bj(q) by the product sin [Tr(qj — %)) sin (jr(qj — /%)), and gives a function 
0{h). Taking the error into account, we obtain \\gj IhWipq^) — 0(h). 

In the case of g^-, i ^ j, we apply sm.{irhDi) to bj(q): the resulting function takes values 
±1 on its support Si fl Sj, so that \\g^ i_H"||L 2 (K™) = 0{h). 

We have now proved that \\F W w^||l 2 (r™) = 0{h) for any F E J. The procedure can be 
iterated to any finite product of functions F« e yielding the estimate (|4.18|) . 

The proof is completed by the periodization argument as in the proof of Proposition 14.81 
The only slight difference lies in the fact that the analogues of the functions u^\q) of 

()4.17|) may now have discontinuities near 3D, so that ||Pt 2 ™ — ^u h ||i2(Rn) = 0{h) 

instead of 0{h°°). □ 

4.3. Quantum relations. Suppose that A C T 2n x T 2n is a Lagrangian submanifold. The 
basic example is given by the twisted graph T' K of a symplectic diffeomorphism k on T 2n 
(see Section 13 .2 .ip : 

K = { (<i',q;p',-p) ■ (q>p) = k (<?,p), (q,p) e T 2n } . 

As we noticed in that section, the choice of change of sign depends on the choice of the 
splitting of variables (q,p), which is itself related with the choice of a polarization in the 
quantization a \— > Op h (a) 24, §25.2]. This somewhat cumbersome convention is explained 
as follows. 

Any state v G Ti.^ is naturally identified to a linear form f v e (Tify* through f v (w) = 
(v, w). In our notations 5 , this scalar product is antilinear in the first component. To make 
the identification linear, we choose instead 

(4.25) ve?%=> /„(•) = (v,*), 

where states v are written as vectors in the basis f)4.3|) . 

Let CiTC^) ~ Hh®(Hh)* be the space of linear operators on TC%. The linear identification 
(EPS} of HI with (HI)* gives the identification, 

(4.26) C(Hl)~H 2 t?, through (u ® v)(w) = u (v, w) , u,v,weH r y\. 

5 This is the physicists' convention. 
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We observe that the norm on Tif^ 1 is the same as the Hilbert-Schmidt norm on C(TC^): 

(4.27) \\T\\ H ^ = (ti H n(T*T))K 
It is related to the operator norm on C(7i^) as follows: 

(4.28) — \\T\\nl n — ^ 11^11 ■ 

In particular, unitary operators have Hilbert-Schmidt norm N n ' 2 = (2vr/i)- n / 2 . 

The identification ()4.26|) dictates the way an operator of the type A\ ® A 2 (with Aj G 
C(H^)) acts on u ® v G 7i 2 ™ ~ C(H]l). Indeed, if we take any w G H^, we have 

[(Ai <g> A 2 ){u g) «)](«;) = [Aim ® A 2 u](w) 

= Aim (A 2 u,«j) 

= Aim (u, A 2 w) 

= [(Aiu <S>v)o A' 2 ](w) . 

Here A' 2 is the transposed of the operator A 2 , written as a matrix in the basis ()4.3j) . In 
the case = Op h (a), A 2 = Op h (6) for some real functions a,b £ C°°(T 2n ), one checks 
that A' 2 = Op h (b), with the same twisted function as in the proof of Lemma [4.61 b(q,p) = 
b(q, —p). By linearity, for any Ch G TC^ 1 ~ C(H^), we have 

(4.29) Op h (a ®b)C h = Op h (a) o C h o Op h (b) . 

The sign change in the tilting r ~> V parallels the transformation a(p') b(p) ~» o(p') 

We are now in position to quantize a symplectic map, more generally a symplectic relation 
T as defined in Section l3~2l 

Definition 4.11. A semiclassical sequence U — {Uh G 7i^ n } h _> satisfying 

(4.30) \\U h \\ n 2n = 0(h K ) , for some fixed K G R, 

is a quantum relation associated with the symplectic relation r if U is a Lagrangian state 
in I(r'), in the sense of Definition 14.71 

Explicitly, for any M G N and any sequence of functions 

9 3 e C°°(T 2n x T 2n ) , ^|iv = 0, 1<J<M, 

we must have 

(4.31) || Op^Af) O • • • O OpftGft) C^H^ = 0(/l M ) ||C/J| W 2n . 

The assumption that [7^ is tempered in the sense of ()4.3U|) (which also implies tempered- 
ness in the operator norm) is necessary to assure that composing Uh with residual (0(h°°)) 
terms produces residual terms. That is a standard assumption in C°° semiclassical calculi 
- see PQ EI], and will be used in the proof of Prop J4. 121 The quantum weighted relations 
defined in § 14.41 will naturally be tempered, having norms ||?7/ 1 || W 2n = 0{hr n l 2 ). 
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If a function g G C°°(T 4n ) vanishes on F, then the function g defined as g(q',p'', q,p) = 
g(q' ',p'; q, —p) vanishes on V. The condition gr 3 -|p = can thus be written gj\r = 0. 

We also note that (j4.31j) entails a version of Egorov's theorem. If fi, /r G C°°(T 2n ) 
satisfy 

(p'p)er =► f L (pr) = f R (p) , 

then we have 

(4.32) || 0p,(/ L )f4 - U h O Vh (f R )\\ Hl n = 0(h) \\U h \\ Kn . 

Indeed, the function / = f fi <8> 1 — 1 <E> fn vanishes on T, so that / vanishes on F. We 
then simply apply the definition (|4.<Uj) with g± = f and use ()4.29j) . When T is a graph of 
a symplectic transformation, fn is the pullback of fi, and we get a statement similar with 
the standard Egorov's theorem. 

Remark 4.2. Following Section [Ol in the case when F is a Lagrangian with boundaries 
projecting on a hypercube, it is useful to include in the definition sequences U = { Uh } in 
the (larger) space I(F' U A s ; ng ); the quantum baker's relation we define in next section will 
belong to such an enlarged space. 

Through the identification (|4.26jl . Uh is an operator on Ti%. We now show that this 
operator "classically transports" the microsupport of a sequence w = { Wh E Ti.^ }. 

Proposition 4.12. Take U = {U h G H 2 h n ~ C{Hl) } a quantum relation U G 1(F). Then 
for any sequence w = {wh G Ti.^}, \\wh\\ ^ 1, the microsupport of the image sequence 
U(w) = { Uhiwh) } satisfies: 

WF h (U(w)) C r(WF h (u>)) = { p' G T 2 " : 3 p G WF h (w) , (p',p) G T } . 

Proof. Assume that p' £ T(WF h (w)), which means that r _1 (p' ) n WF h (w) = 0. Then 
there exists a function / G C°°(T 2n ) with / = 1 near Pq but with supp(/) sufficiently small 
so that T _1 (supp(/)) d CWF^w). Consequently, there exists a function g G C°°(T 2n ) 
with g = 1 near WF^(w) but g = on T _1 (supp(/)). The function / ® g G C°°(T 4n ) then 
automatically vanishes on F. 

Our aim is to show that p' G" WF h(U(w)). For this, we introduce one further function 
a G C°°(T 2n ) such that a(p' ) > and / = 1 on supp(a). As in the proof of (|4.9|) we see 
that for any M G N, Op h (a) Op ft (/) M = Op fc (a) + C(/i°°). Hence 

|| Op,(a)£W|| = || Op h (a) Op,(/) M U h w h \\ + 0(ft°°) 

< ||Op,(a)0 P?i (/) M ^Op,(g) M u; /t || 

+ || Op A (o) Op,(/) M % (1 - Op h (g) M ) w h \\ + 0(h°°) 

To bound the second term on the right hand side, we notice that the function (1 — g M ) 
vanishes near WF^jo), so from Lemma l4~B1 we get ||(1 — Op h (g) M ) Wh\\ = 0(h°°); from the 
temperedness of Uh, the second term is thus residual. 
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The first term on the right hand side is estimated using the identity 

Ph (f) M U h Ph (g) M = Op & (/ ® ~g) M U h . 

Because f<S>g vanishes on T', the Hilbert-Schmidt norm of that operator is 0(h M+K ), where 
K comes from the temperedness of Uh, f!4.30|) . Using ()4.28|) . we thus get || Op h (a)UhWh\\ = 
0{h M+K ) for an arbitrary M G N, which shows that p' Q £ WF h (U(w)). □ 



4.4. Quantized weighted relations. In Section 13.41 we equipped a symplectic relation 
T with a measure, or weight fi. In order to associate to the weighted relation (T, fi) a 
sequence of operators Uh G 7"^™, we need to elaborate on Definition 14.111 thereby defining 
a subfamily I{T',fi) C J(T'). 

In the standard microlocal context |2H Section 25.1], a Lagrangian state ip G /(A) has 
a well defined amplitude, or symbol, which is a section of the Maslov half density bundle 
over the Lagrangian submanifold — see [23 Theorem 25.1.9]. The local aspects of this 
procedure have recently been adapted to the semiclassical case PQ, and a similar approach 
can be used in the case of T 4n . 

Although one could characterize the operators quantizing (T, /z) in terms of their symbols 
(grossly speaking, the absolute square of the symbol should equal the weight jj), we won't 
do it here, in order to avoid technical issues involved in the description of the symbol map. 
Instead, in the definition below we use bilinear expressions in Uh, which allows us to avoid 
introducing symbols. 



Definition 4.13. Let (T,/i) be a weighted piecewise smooth relation as defined in 
and let U G I(T' U A sing ), in the sense of Definition 14.111 and Remark 14.21 For any \ a G 

C°°(T 2n ; [0, 1]), a = L, R, we define 

U XLXR = °Ph(XL) U h OphiXR) ■ 

We say that U quantizes the weighted relation (T, //) if for all xl, Xr with sufficiently small 
supports satisfying supp(x L <g> xr) H A' sing = 0, 

u ^ u XLXR u* XLXR = 0Mgl LXR ) + o(h) 

' U x,Xn U x,xn = Oph(9 X R LXR ) + 0{h) , 

where g^ LXR are the functions given in (|3.9|) . and the remainder is 0(h) in the operator 
norm on C(Ti.^). We then write 

u = {u h }ei(r'uA sing ,p). 

The conditions on the smallness of supports of Xa guarantee that the operators appearing 
on the left in ltP3|l are of the form Op h (/), / G C°°(T 2n ). That follows from the fact that 
T is locally a graph — see < J3.4I 
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If T is the graph of a symplectic diffeomorphism k and p = ■K* L {u n / n\) = ir* R (u) n / n\) , then 
Uh is unitary to leading order: 

U* h U h = I + C h , U h U* = I + D h , \\C h \\ c{Hl) = 0(h), ||D ft |U ( «n) = 0(h) . 

For h small, (/ + Ch)~~ z , (I + Dh)^ 1 exist, therefore a possibility to make the quantization 
strictly unitary is to replace Uh by UfJJ + C^) _ 2 or (J + D h )~^Uh- 

The condition (|4.33j) can be interpreted as follows. Suppose that if) G ||^|| = 1, is 
micro localized at a single "regular" point p : 

WF h (iP) = { Po }cT 2n \ 7r R (A' aing ), 
and r(p ) = U/ =1 Pj-, = /%(p )- Then, 



3=1 

ll^-f = P,(p ) + 0(h) , WF h (^) C { p' 3 } . 

From Lemma 14. 5( if Pj(po) ^ then WF/,(^-) = { }• A similar statement holds for U£. 

Indeed, if for each j = 0, • • • , J we take Xj G C°°(T 2n ; [0, 1]) supported in a small neigh- 
bourhood of p , resp. Pj, and equal to 1 near that point, (j3.10p shows that g^ xo (p, 
Pj(po) for j = 1, . . . , J. On the other hand, Proposition 14. V2\ gives 



0J 



3 ' 

J 



U h if) = U h Op,( X o)^ + o(h°°) = Y, u XjX0 V + 0(0 , 
l 4 -^ 4 ) j=l 



wF,(f/ XjX0 ^)c{p;} 



If we take if>j = f U XiXo ip then 



J w XjXo 

2 - (^xo^W, V>> = (0 P/ K JXO )^> + = Pi(Pd) + 0(h) . 

Example. We now consider a special case of quantum relations Uh, of the form 

(4.35) (Qj\Uh\Q k ) = N- n / 2 a(Q h Q k ) exp (2mNS(Q J ,Q k )) , 

where a, S 6 C 00 (T 2ri xT 2n ) and the generating function S^g', g) satisfies the nondegeneracy 
condition det(<9 2 ,^) ^ near the support of a(q',q). Using Definition 14 . 1 II we see that Uh 
is associated to the graph T s of the symplectic transformation (q, —d q S) i— > (q',d q 'S). To 
be more precise, 

(4.36) U h e I(T' s ,/i s ) , for p 5 = f |a(g', g)| 2 dg' dg , 
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where we used the coordinates (q', q) on T s . Projecting this measure on the left and right 
tori, we get: 

(4.37) 

q:p'=d q ,S(q',q) 

The above sums are always finite. This example will be used to analyze the quantum 
baker's relations studied in the next sections. 



4.5. Quantized baker's relation. We explicitly construct quantum relations e £{7~L\) 
associated with the "open baker's maps" described in Section I3~3l For simplicity, we will 
assume that the coefficients Dj and £j are integers. Besides, we will only consider the 
subsequence of Planck's constants of the form h = (2ttN)~ 1 such that N/D\ = Mi e N 
and N/D 2 = M 2 G N (that is, iV is a multiple of lcm(£>i, D 2 )). 

Restricting ourselves to this subsequence, we define the quantization of the baker's re- 
lation (|3.4|) as the following operators (written as iV x N matrices in the bases (|4.3Jl ): 



(4.38) 



D def ^ 
&h — J~K 



N 



(0 











o\ 












































2,h- 



The numbers of columns in successive blocks are respectively given by 

£ 1 M 1 , M 1 , £ 2 M 2 - {t x + \)M X , M 2 , (D 2 -£ 2 - l)M 2 , 

and Tm is the discrete Fourier transform given in (|4.5|) . These matrices obviously generalize 
the unitary matrices associated with the closed baker's map [2]. 

We now check that the matrices ()4.38|) satisfy the Definition 14.131 if we select the appro- 
priate Lagrangian surface on T 4 , namely by adjoining a singularity set A sing to the twisted 
graph B' (see Remark l4.2j) . and equip B with the weight fi described in (j3.11|) . By linearity, 
we can separately consider the two blocks Bj )h . Let us study the left block Bi h . Since the 
classical relation Bi is generated by the function Si(q,p') of (|3.5|) . it is natural to express 
the operator B^ in the mixed representation (p',q), that is by a matrix from the basis 
{ \Qj) } to the basis { \Pk) }. Since the change of basis matrix, (\Pk) (Qj\)j,k=o,...,N-i, equals 
JF/v, the operator defined as the matrix 



((Qk\Ai,h\Qj))j,k=o,...,N-i — {(Pk\Bi t h\Qj))j,k=o,...,N-i — Fn ° Bi t h 

is given by the Fourier block J r Ml at the same position as in (|4.38|) . and zeros everywhere 
else. 
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The following lemma reduces finding the (weighted) Lagrangian relation associated to 
B\ h to finding the (weighted) Lagrangian associated to A\^. We denote by F the following 
transformation of T 2n : F(q,p) = (p, —q). It means, we rotate by —n/2 around the origin 
in each plane (qi,Pi)- We denote by Fl the transformation of T 4n acting through F on the 
left coordinates (q',p') and leaving the right coordinates unchanged. 

Lemma 4.14. Suppose that Uh G £(7^) — 7"^ n and that Vh = Fn Uh- Then, for any 
(possibly singular) Lagrangian C G T An , 

U h g 1(C) ^V h e I{V) , 

where 

V = F L (C) , equivalents V = F L (C) = {{jp' , -q'; q,p) : (q',p'; q,p) G C} . 
Furthermore, 

U h El(C',fx) ^V h eI(V',v), with v = F l * I m. 

Proof. The transformation C — > T> results from a general composition formula which can be 
proved by mimicking the semiclassical proof in £Q . Here it follows from the the covariance 
properties of Weyl quantization with respect to the Fourier transform: for any a G C°°(T 2 "), 



(4.39) T h x Op h (a) ° T h = Vh (a o F) . 
As a result, for any / G C°°(T 4n ), 

Op h (/)(^ A o U h ) = T h o Op h (/ o F L ){U h ) . 

This identity proves the first assertion. 

Using ()4.39|) . we notice that for any xl, Xr G C°°(T 2n ; [0, 1]), the cutoff propagator 
V XlXr satisfies 

VLx* v xlx r = U* XLoFxR U xloFxr = Op h (g x R L ° FxR ) + O(h) , 

V XLXR V; LXR = T h U XloFxr U* XL0FXR Tl = Op h (9 X L L ° FXR o F- 1 ) + O(h) . 

Using the pushforward of functions Fl* f = f ° F^ 1 and the fact that itroFl = ttr, we get 
9r L ° Fxr = ^r*X^*l( f l*Xl) XrXr y) = ^r*(Kxl KrXr F l ^) 
g x L oF XR Q F -i = k l *F l *{-kI{F£ X l) itrXr V) = *L*«XL TTrXR F L *fi) . 

This proves that Vh is associated with the weight v = Fl* fionV. □ 

Let us now describe the weighted Lagrangian associated with the operator A^. The ker- 
nel of that operator vanishes outside the square H — Ii xix, where Ji = [£i/D 1 , l\ + 1/-Di], 
and on H it takes the values 

(4.40) (Q k \A lth \Qj) = {P k \B^ h \Q 3 ) = J^ l H (Qk,Qj) exp ( - 2iirN S 1 (Q k ,Q j )) . 
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The operator has the same form as in ([4.35)1 . with the (obviously nondegenerate) 
generating function S = —Si and symbol a(q',q) = \fD\ ^h{q', q)- If we forget (for a 
moment) the discontinuities of the symbol, we find that A\^ is associated with the graph 

r Sl = { (q\ ~(D iq - 4); q, (D iq ' - l x )) : q, q' G h } , 
equipped with the weight 

fx Sl = Di l H {q\q) dq dq. 
From Lemma f4. 14) the operator Bih = JF^ o A^h is associated with the graph 

F^(T Sl ) = { {{D x q - 4), q'; q, {D 1 q > - £,)) : q,q' G h} = B, 
and the weight 

A*i = f l*VSi = Di l H (p, q) dp' dq , 

which can be expressed as 

/ii = t h (q) dq dp , 7r L * /ii = l h (p') dq' dp' . 
It represents the half part of the weight (|3.11|) . 

Let us now take the discontinuities of a(q', q) into account. Since they occur at the 
boundary of the square H, they have the same consequences as in Lemma \4. 91 Namely, we 
must add to the Lagrangian T' s a "singular" Lagrangian, which is the union of 4 pieces, 
each piece sitting above a side of H. This Lagrangian should then be rotated through F^ 1 
as well. 

For instance, the side {q' G I\, q = t\fD{\ leads (after rotation) to the singular La- 
grangian 

Asmg.i = | (q = 0, q = jj-; p' ,p) : p' G h ,p G 1 1 , 

which contains the corresponding side of dB[. Similar Lagrangians A singi j, i = 2,3,4, 
contain the other sides of dB[. 

The same analysis applies to B 2) h and hence we have proved the 

Proposition 4.15. The sequence of matrices { Bh } given in ()4.38|) quantizes the clas- 
sical baker's relation B = B\ U B% of ()3.4|) . in the sense of Definitions \^~TJ[ \4-l§ and 

Remark 

B h el(^B'u\jA 

where the weight /i is given by (|3.11|) . 

This quantization of the baker's relation is very close to the "quantum horseshoe" defined 
by Saraceno-Vallejos in The operator B^ is contracting, and its eigenstates can be seen 
as "metastable states", "decaying states" or "resonances". This contraction mirrors the 
decay of a classical probability density evolved through the open map B (due to the "escape" 
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of particles to infinity). This classical decay can be analyzed in terms of a "conditionally 
invariant measure" on T 2 jH], which decays according to the classical decay rate 7^ = 



4.6. Numerical check of the Weyl law for the baker's relation. We have numerically 
computed the spectra of the quantum baker relations for various symmetric and nonsym- 
metric baker's relations. Results for the symmetric "3-baker" (D — 3, 1 — 0) were presented 
in (see also Table [I]), some for the "5-baker" (D = 5,£ = 1) were given in [3U|, while a 
nonsymmetric map (D\ = 32, D2 = 3/2) was studied in [HH]. In the symmetric cases, the 
trapped set is a pure Cantor set of dimension 2d = 2j^jy, so that for any 1 > r > 0, the 
number of resonances in the annulus { |A| > r } is expected to scale as 



in the limit N — ► 00. Our numerics for both maps shows that this scaling is roughly 
satisfied along any sequence N — > 00; much better fits are obtained for N taken along 
geometric sequences of the type iV = N Q D k , with iV fixed and k — > 00 (as in Table Q), 
which lead us to the following weaker conjecture for the symmetric maps: 



#{A G Spec(£ fc ) : |A| > r} ~ C(N Q , r) , (27ih)- 1 = N = N Q D k , k -> 00 . 



Here, the "profile function" C(N a ,r) may (slightly) depend on the "root" of the geometric 
sequence. The special role played by geometric sequences is probably due to the strong 
relationship between the symmetric D-baker and the D-nary decomposition. 

On the opposite, for the nonsymmetric map the fractal Weyl law seems accurate for an 
"arbitrary" sequence iV — > 00 which was also the case for the nonlinear map studied 



Let us explicitly compute the matrix elements of the two vertical blocks Bi ^, B 2} h i n 
(|4.38|) . for the symmetric 3-baker. Both are matrices iV x N/3, which we index by < k < 
N- 1, < I < N/3 - 1: 



log(L>r 1 + J D 2 - 1 ). 



#{A e Spec(B h ) : |A| > r} ~ C(r) N^ 



by 02| . 



5. A TOY MODEL 



(5.1) 




luI)/N if k ^ 31, 
ifk = 31, ' 



(B2,h)kl — ^3 (Bi t h)kl , 



where u)% = e 



2ivr/3 



The largest matrix elements are near the "tilted diagonals" k ~ 31, and decay as 
1/1 A; — 3l\ away from them (see the Figure 6 in |1Q|). Being unable to rigorously analyze 
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the spectrum of Bh, we replace this matrix by the following simplified model: 

Bh = B N = [Bi t h, 0, B 2> h] , 
1/V3 if/=|*/3j , 5 > , M 



(5.2) 



kl 







if I ^ [k/3\ 

where [^J denotes the integer part of x. For iV = 9, this gives 



B2,h)kl — <^3 (B\ : h)kl , 



(5.3) 



Bn =9 — 



1 

v 7 ^ 



/ 1 

















1 





\ 


1 

























1 




























1 

















1 








1 

















w§ 








1 




























1 

















1 








1 

















^3 / 


\o 





1 


















, ^3 



3 27ri/3 



The model has been obtained "by hand", by replacing "lower order" terms in the ma- 
trix Bh by 0, keeping only nonzero elements on the "tilted diagonals", and replacing 

(1 _ e 2 ^( ±1 )/ 3 )/(AT(l - e 2 M±i)/^)) by 1. 

The new matrix Bh retains some qualitative features of Bh but there is no immediate 
connection between their spectra: the "lower order" terms are not small enough for that, 
and Bh cannot be considered as a "small perturbation" of Bh,- 

The simplicity of the matrices Bh will allow us to prove (in the case iV = 3 k , k e N) 
the fractal Weyl law which we could numerically observe for Bh (see Section f5.2|) . It is 
interesting to notice that the simplified operator Bh is in fact not associated with the same 
classical relation as B^: 



Proposition 5.1. In the notations of Section \4 . 4 the quantum relation {Bh} is associated 
with the weighted relation (B,p) given by (I3.12j) and f)3.13|) : 



■A-sing 



j=0 



B h G I(B' U A sing , jl) , where 

{ (Y = 0, q = j/3 ;p',p), p', p E I } . 



A-singJ 



Proof. In place of Bh we will consider Ah = o Bh, and apply Lemma 14.141 From the 
structure of Bh, the operator Ah can obviously be split into + ^2,h- We will analyze 
the first component in detail, the analysis for the second one being similar. The matrix 
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(Qk\^i,h\Qj) is nonzero in the vertical strip I x x I, with Ii = [0, 1/3): 



(Q k \AiM) = (]Te- 2 ^) expi-QiirNQuQj) . 

V3iV <=Q 

Like A\ t h (see ^4.5|) . this operator is of the form ()4.35|) . with generating function S(q', q) = 
—Si(q',q) = Sq'q and discontinuous symbol 

e- 2 ^' sin(37rg') 

a g ,g) = l 7l (g) — -= — — — . 

y/3 sm(7rg / ) 

Forgetting about discontinuities, A\^ is therefore associated with the graph 
r Sl = {(q',p' = -3q;q,p = -3q'), : q' G I, q G I x } , 

and the weight 

i / / \:2 i i j n / \ sin 2 (37rg / ) , 
= a(g , g) dq dq = l h (g) . 2 — - dg dg . 

dsm (vrg'j 

After applying the transformation of Lemma 14.141 this leads to the graph 

2 

^(fsJ = {(</ = 3g, g;p', p = 3p') : g g 7 1s p' g I } = |J B y , 

3=0 

and the weight 

-it / \ sin 2 (37rp') , . , 
i^ 1 ft Sl = t h (q) , \ P dp' dq . 

dsm [Tip) 

Through the change of variable (g,p ; ) i— > (g,p), we see that this is the weight ()3.13j) on the 
component B\. 

The discontinuities of a(q f , q) only occur along the two segments {(q' G I, g = 0)}, 
{(g' G I, g = 1/3)}: they generate the singular Lagrangian 

= j = 0, q = J -; p' G I, p G I) j , j = 0, 1 , 

which transforms under F^ 1 into the components A sing , A sing )1 . 

Similarly, the second part of the matrix, I?2,/i, is associated to the twisted graph B' 2 with 
weight /2,g 2 and the two singular components A singj2 , A sing0 . □ 

As explained in Section 13.41 the graph B can be obtained by adjoining to each point 
(p'; p) G B the points (p' + (0, 1/3); p) and (p 1 + (0, 2/3); p). This "aliasing" is due to the 
diffraction created by the sharp cutoff in the matrix Bh, as opposed to the "smooth" decay 
of coefficients in Bh- A similar aliasing was observed in |HH] for the graph associated with 
the unitary matrices A 2 k defined in (15. 9)) : instead of quantizing the standard 2-baker (|3.2|) . 
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they are associated with a multivalued map obtained from it by aliasing. This observation 
was obtained using the propagation of coherent states. 

Both B and B share the same forward trapped set r_ = T_ = C x I (see fcj.3|) . but the 
backwards trapped set of B is easily shown to be T + = T 2 , which drastically differs from 
T + . This asymmetry between T_ and r + reflects the fact that, unlike B, the relation B is 
not time reversal symmetric. 

The fact that Bh is not associated with the relation B should not bother us too much 
though. In the next section, we will give a more "formal" construction of the matrix B^, 
in the case where N is a power of 3 (this construction will also hold for any symmetric 
.D-baker, for N a power of D). We will show that this matrix naturally appears through a 
"nonstandard" (Walsh) quantization of the open 3-baker relation B. 

5.1. Walsh quantization of the baker's relation. The Walsh model of harmonic anal- 
ysis has been originally devoted to fast signal processing It has been used recently in 
mathematics to obtain simpler (and provable) versions of statements of the usual harmonic 
analysis — see jSSI for an application in scattering theory and for pointers to the recent 
literature. The major advantage of Walsh harmonic analysis is the possibility to completely 
localize a wavepacket both in position and momentum: for our problem, this has the effect 
of avoiding diffraction problems due to the discontinuities of the map, which spoil the usual 
semiclassics j3Hl- Closer to our context, Meenakshisundaram and Lakshminarayan recently 
used the Hadamard Fourier transform (which is related with the Walsh transform we give 
below) to analyze the multifractal structure of some eigenstates of the (unitary) quantum 
2-baker B h (23). 

5.1.1. The quantum torus as a system of quantum Dits. We first fix the coefficient DsN 
(D > 2) of the symmetric baker's relation (J3.8|) . and will consider in this section only the 
inverse Planck's constants of the form N = D k for some k G N. In this case, integers 
j G 1> D k = {0, . . . , D k — 1} are in one-to-one correspondence with the words e = eie2 •••£*; 
made of symbols (or "Dits") q G 1>d'- 



The natural order for j G Z D fe corresponds to the lexicographic order for the symbolic words 
{e G {TjdY}- This way, each position eigenstate \Qj) of the basis ()4.3|) can be associated 
with the unique symbolic sequence 6162 • ■ • which gives its Dnary expansion 



Let us denote the canonical basis of C D by {eo, e±, . . . , &d-i\- Then, each \Qj) can be 
written as 



a.- 



(5.4) 




(5.5) 



Qj = — = 0- t x e 2 ■■■e k . 



(5.6) 



\Qj) = e ei <8> e €2 <g> • • • <g> e, 
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Following [48 , we denote each \Qj) by |e) = |eie 2 • • • e^) to emphasize the above tensor 
product decomposition. This way, the quantum space TL\ is naturally identified with the 
tensor product of k spaces C D : 

n\ = (c D )i® (c D ) 2 ®---® (c D ) k . 

In the quantum computating framework, each space (C D )g is interpreted as a "quantum 
.Dit", or " quDit", and the basis { |e) } is called the computational basis [35J. Viewed in 
our toral phase space, the qu.Dit (<C D )e is associated with the scale D~ e in the position 
variable, so (C D )\ is called the "most significant quDit". 



5.1.2. Walsh Fourier transform. The discrete Fourier transform of (J4.5)) (with n — 1, N — 
D k ) is the Fourier transform (in the sense of abstract harmonic analysis) on the group Z D fc. 

More explicitly, each row of J- D k corresponds to the character f \— > exp ^ — 2iirjj'/D k ^j of 

Z D fe. Using (J5.4)) . the matrix elements can be factorized: 



D e 



(5.7) (F Dk ) 3f = D- k ' 2 exp ( - 2i7r^) = D^ 2 JJexp ( - 2i 

Notice that each e m (jj') can be easily expressed in terms of the symbols of j and j': 



The Walsh Fourier transform is the Fourier transform on the group (Z^) fc . It can be defined 
by keeping only the first factor on the right hand side of (|5.7|) : one obtains the matrix 

(5.8) (W k ) jf = D- k ' 2 exp ( - 2z7T^^) = J] D' 1 ' 2 ^H+i-tW) ^ ^ = e& /fl 

£=1 

Using the identification 7i\ ~ (C D )® fc , this definition can be recast as follows. 

Lemma 5.2. The Walsh Fourier transform W k acts simply on tensor product states: 

Wk(vi (g) • • • ® v k) — FoVk ® ■ • ■ ® FdV\ , ve e C D , £ = l,...,k. 

Here Tb = W\ is the discrete Fourier transform on C D . As a result, W k is a unitary 
tranformation on Ti^. 

The proof consists in a straightforward algebraic check. 

As opposed to the discrete Fourier transform, the Walsh Fourier transform does not 
entangle the different qu-Dits: a tensor product state is sent to another tensor product 
state. 
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Example. To illustrate this simple lemma we take D = 2, and consider the following 
2 k x 2 k matrix, k > 1: 



(5.9) A 2 k = [v4 0i2 fc,v4 li2fc ] , {A j>2 k] 



For instance when k = 2 we get 



0<n<2 fc -l,0<m<2 fe - 1 -l 



(-l) jn /V2, m = [n/2\ 
0, m ^ [n/2\ . 



A 2 2 = —= 



1 

v 7 ^ 



A 
1 
1 

V° 1 



i 

-i 






o\ 



1 

-v 



This sequence of matrices has been obtained as the "extreme" possibility among a family of 
different quantizations of the (closed) 2-baker's map [IE! 6 , and its semiclassical properties 
were further studied in (SSj. In a different context, this (unitary) matrix belongs to the 
family of transfer matrices associated with the de Bruijn graph with 2 k vertices [53] . 

The transformation A 2 k acts as follows on tensor product states: 

Vi® ■ ■ ■ ®Vk i — > v 2 <g> • • • <8> v k <8> ■ 

This implies that this matrix can be easily expressed in terms of the Walsh Fourier transform 
(for D = 2): 



(5.10) 



A, 



W, 



k-1 







w, 



k-1 



where the 2x2 block structure corresponds to the most significant (leftmost) qubit. This 
expression exactly parallels the one defining the Balazs-Voros (unitary) quantum baker j2|. 
Compared to this "usual" quantum baker, A 2 k is thus obtained by replacing the discrete 
Fourier matrices T 2 k, T 2 k~\ by their Walsh analogues Wu-i- 

The matrix A 2 k is unitary; as we will see in the next section, our toy model for the 
quantum open 3-baker (see Eq. (|5.2j) ) is its subunitary analogue. 



5.2. Resonances for the Walsh quantization of the open baker relation. In this 
section we set D = 3, and concentrate on the symmetric 3-baker (|1.3|) . By analogy with 
the example in last section, we modify the quantization (j4.38l5.1|) . in the case iV = 3 k , by 
replacing the discrete Fourier matrices by their Walsh analogues. The resulting operator 
exactly coincides with the toy model ()5.2j) introduced in the beginning of this section: 



'We thank M. Saraceno for pointing out this interpretation to us. 
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Lemma 5.3. In the case N = 3 k , the matrix B h defined in (|5.2|) can be rewritten in terms 
of the Walsh Fourier transforms as follows: 

B h 













w* k | 













V o 








We omit the simple algebraic proof. If we define the "truncated" inverse Fourier matrix 

'10 1 

(5.11) Tt = ^=|1 ul 



- 1 h „ 
V 3 \1 cj 3 



the toy model acts as follows on tensor product states: 
(5.12) B h (v! ® • • • g> v k ) = v 2 <8> v 3 <g> • • • ® JJui . 

This form is particularly nice to compute the spectrum of 5^. We start by computing the 
spectrum of its power (Bh) k , which is enough to obtain the radial distribution of resonances 
(that is, the distribution of resonance widths). 

Proposition 5.4. Let X±, |A_| < |A + | ; be the eigenvalues of the matrix 

The non-zero eigenvalues of (Bh) k (for N = (27rh)~ 1 = 3 k ) are given by X k [ T p X p _, < 
p < k, each occurring with multiplicity . From this we get the radial distribution of the 

eigenvalues of B h (counted with multiplicities): 

Vr G [0, 1], 1 # { A G Spec(B fc ) : |A| > r } ^ C(r) , 

(5.13) 

C(r) = 



r < | det Q 3 | 2 
r > I det I 5 



Hence the nontrivial resonances accumulate near the circle of radius r (B) = | detfi3|2. 

This proposition gives Theorem ^ where B is the baker's relation described in Propo- 
sition 15.11 Bh the matrices ()5.2|) . and Planck's constants are taken along the sequence 
{h k = (2tt x 3 k )-\ k G N}. 

Proof. From the expression ([5.12)1 . we see that 

{B h ) k {v x <g> • • • ® v k ) = F£v x ® • • • ® ^ 3 *ffc • 

That means that {Bh) k = (-7"|)® fc , so one eigenbasis is obtained by taking the tensor 
products of eigenstates of and the eigenvalues of (B h ) k are the corresponding products 
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of eigenvalues of T^. The nonzero eigenvalues A + , A_ of Tt, are the eigenvalues of f2 3 , so the 
first part of the proposition follows. To prove the second part, notice that each eigenvalue 
y^-vyv_ Q f correS p 0nc l s to an eigenvalue (possibly in the generalized sense) of modulus 
| A i_-p/* a p/*| o{ Q h Therefore, we are able to count eigenvalues of Bh (with multiplicities) 
in a given annulus. 

Let H(t) denote the Heaviside function, H{t) = for t < 0, and H(t) = 1 otherwise. 
Then, for any < r < 1, 

# { A e SpecCB,) : |A|>r} = ^iJ(|A + | 1 -^|A_r-r)( A; ) 

p=0 

E, , n x A \ loerfl A_A + |i/r) 

p J { - p/k+1/2 + P) U' " = lo g (|A + |/|A.|) • 

Using Stirling's formula, one easily gets in the limit k — > oo: 

J_ £ H ( _ p/t + 1/2 + p) Q „ y| £ e — * ix _ HW . 

This expression shows that the distribution of resonances is semiclassically dominated by 
the degrees |p — k/2\ = 0(k 1 ^ 2 ), and proves the second part of the proposition. □ 

The explicit eigenvalues are A± = ^^jp ± ^/ n-^^ w ith approximate values 
A+ « 0.8390 + ^0.0942, |A+| w 0.8443, A_ « -0.5504 + i4058, |A_| « 0.6838. 



The geometric mean of their moduli is r (B h ) = lA-A+j 1 / 2 = a/| det f^] = 3 1 ^ 4 . 

We need to analyze the spectrum of 5^ more precisely to show that the distribution of 
resonances is asymptotically uniform with respect to the angular variable. 

Proposition 5.5. Let h = (27r3 fc ) _1 . As a set, the nontrivial spectrum of Bh is given by 

{A+} U {A_} U |J {u\^ plk \ p l k : 1 < v < k - 1} . 

UJ k = l 

For each p 7^ 0, k, the k eigenvalues asymptotically have the same degeneracy | , which 
shows that their distribution is uniform in the angular variable. Therefore, for any contin- 
uous function f e C(D(0, 1)) we have (counting multiplicities in the LHS): 

0^AeSpec(S h ) 
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Proof. To classify the nontrivial spectrum of Bh, we will use the eigenvectors v± of F£ 
associated with the eigenvalues A±. Call 

{ 1 = 77i?72 • • • Vk ■ V£€ {±1 } - ( z 2) k 

the set of binary sequences of length k. The number of symbols % = — in the sequence 
■q is called the degree of 77. The cyclic shift r acts on these sequences as t(t]i • • = 
V2 • ' • VkVi ■ The shift allows us to partition (Z 2 ) fe into periodic orbits, each orbit O = 
{ rj, r(rj), . . . , r £o_1 (77) } being of (primitive) period £q = £ v - Since r k = id, the primitive 
period must divide k. We call deg(O) the common degree of the elements of O and observe 
that 

(5.14) k\£ deg(0). 

To each sequence r) we associate the state \rj) = f v m ®v m ®- ■ -<S>v Vk , which is obviously an 
eigenstate of (Bh) h , with eigenvalue x^-deg^^deg^ ^j^gg states form an independent 
family, which span the nontrivial eigenspaces of Bh. This operator acts very simply on 
these states: 

Vr7G(Z 2 ) fc , B h \rj) = X m \r(rj)). 

Hence, for any orbit O, Bh leaves invariant the ^-dimensional subspace Vq == span { \r\) , 77 G 
To compute the spectrum of Bh\v a we fi rs ^ observe that it is contained in the set of k-th 
roots of ^k-dcg(e>)^deg(e>)^ w j 1 j c j 1 j n v j ew Q f (|5.i4j) i s equal to 

S d = f { 4 o A ^deg(a)/ fcA d_cg(o)/ fc ; . = Q) j £o _ j J 

We claim that Spec (5^1^) = Se> (clearly with no degeneracies). In fact, let VL : Vq — > Vq 
be defined by Vto\T l {r])) = Lu^\r e (r})), for a choice of r\ G (9. This operator is invertible 
on Vo. By a verification on basis elements, 

Bh\v a — Bh I v 5 

and hence if A G Spec (Bh\v a ) then ^ C A G Spec(5/ l |y ) for any j. 

Since O ^ O' => Vq fl V©/ = {0}, enumerating the orbits decomposition of (Z 2 ) fc 
yields the full nontrivial spectrum of Bh- In spite of the large degeneracies, this nontrivial 
spectrum does not contain any Jordan block. 

The degree p = corresponds to the unique orbit O = {rj = + + ■•■ + }, so the eigen- 
value A + is simple. Similarly, the degree p = k leads to the simple eigenvalue A_. 

For any degree 1 < p < k — 1, call g = gcd(k,p). The sequences of degree p will take 
all possible periods £ v = kjl, where £ G N, £\g. We show below that, in the semiclassical 
limit, the huge majority of the sequences of any degree p 7^ 0, k have primitive period k. 
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Lemma 5.6. There exists C > 0, K > s.t, for any k > K and any degree 1 < p < k — 1, 

#{ryG(Z 2 ) fc : degfa) = p , £ v < k } < ^logfc 
#{77 e (Z 2 ) k : deg(T7) = p} ~ k 

Proof. We still use g = gcd(/c,p). 

If g — 1, then all orbits of degree p are of primitive period k. 

If g > 1, there exists i > 1, For any P prime divisor of £, any sequence of primitive 
period £ v = k/£ is also of (nonnecessarily primitive) period k/P. Any sequence of degree p 
and (nonnecessarily primitive) period k/P can be seen as the P repetitions of a sequence of 
k/P bits, among which p/P take the value (— ). Therefore, the number of such sequences 
is exactly (pjp) ■ As a consequence, we have 

, , 15) #{T7€(Z 2 ) fc : deg(rj)=p,£ v <k} ^ ©3 

#{rye(Z 2 ) fc : deg(rj)=p} ~ ^ (*) ' 

P prime, P\g \pJ 

We will now estimate each term in the above sum. From the symmetry ( ) = we 
can assume p < k/2. Expanding the coefficient ( ) into 

'k\ k(k-l)---(k-p+l) 



j) J p{p — 1) • • • 1 

we notice that both the numerator and the denominator contain exactly p/P factors which 
are multiples of P. Their ratio gives (^p) , while the ratio of the remaining factors is 

Q (k - 1) ■ ■ ■ (k - P + l)(k - P - 1) ■ ■ ■ (k - p + 1) 
(*/£) " (p-l)...(p-P + l)(p-P-l)--l 



k — p + 1 k 

> > - + 1. 

1 - 2 



Here we used the fact that each factor (k — m) /(p — m) > 1, < m < p — 2, and only kept 
explicit the last factor. The last inequality comes from p < k/2. 

We have obtained a uniform upper bound for each term in the sum of (|5.15|) . By 
standard arguments, there exists K, C > s.t. the number of prime factors of any k > K 
is < C log k, so the number of terms in the sum is < C log k. As a result, (j5.15J) is bounded 
from above by C log k/ (k + 2), which proves the lemma. □ 

This lemma shows that the orbits of period <k have a negligible contribution to the 
asymptotic density of resonances. We can therefore act as if, for any 1 < p < k — 1, each or- 
bit of degree p had period k, leading to the k eigenvalues j lu{ \+ p ^ k \ p J k ) j = 0, . . . , k — 1 1. 
In the semiclassical limit, these k eigenvalues are uniformly distributed on the circle of 
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radius |A+ p ' k X_\ : and each of them has multiplicity f (*)■ This shows that the asymp- 
totic resonance distribution is circular-symmetric, with the radial distribution described in 
Proposition 15.41 □ 

Remark 5.1. Several features of the (nontrivial) spectrum of are very different from 
what one expects for a random sub unitary matrix of size 2 k x 2 fc : the (logarithms of the) 
resonances form a regular lattice, most eigenvalues are highly degenerate, and the radial 
density is a delta function at r Q (B). Actually, the only generic feature seems to be the 
global fractal scaling of the Weyl law, and the uniform angular distribution. 

Remark 5.2. The radial density of resonances is governed by tq(B) = yj det f2 3 |, which 
seems to depend on the subtelties of the quantization. As an example of this fact, in ^H] 
we will consider the open baker's map with D = 4, which we call B, obtained by keeping 
only the second and third strips. It has Lyapounov exponent log 4, and the Cantor set C 
(see EI3.3J) has dimension v = log2/log4 = 1/2. The open map B' obtained by removing 
the first and third strips has the same characteristics. However, if we Walsh-quantize these 
two maps, the spectra of B^ and B' h are very different. These spectra are obtained from 
the eigenvalues of different 2x2 blocks of the inverse Fourier matrix T\. The first map 
leads to the block 




with two nonzero eigenvalues A± of different moduli, so the spectrum of B^ will satisfy the 
fractal Weyl law, and be concentrated around the circle of radius 

r (B) = vTdetQlj = 2" 3/4 . 

In an opposite way, the second map leads to the singular block 

The nontrivial spectrum of B' h then reduce to a simple eigenvalue A + = 1. In that case, 
the Weyl law is singular, corresponding to the profile function C(r) = 0. This qualitative 
difference between both spectra cannot be explained from the classical maps, but is due to 
the phases in the matrix elements of the quantum maps. 

6. Conductance in the Walsh model 

6.1. Quantum transport. In this section, we consider open baker's relations for which 
the "opening" consists in two disjoint intervals, which are supposed to represent two "leads" 
connecting a quantum dot to the outside world. We will prove Theorem El in this setting: 
(HHJ) in gOand flQ in §0 
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The baker's relations defined in ^3.31 can all be seen as truncations of invertible maps on 
T 2 . More precisely there exists an invertible baker's map, k : T 2 — > T 2 , such that the graph 
B = B\ U B2 of the open baker's map is 

B = T K n { (q,p) : q G h U J 2 = I, p G 1} . 

For admissible values of N, one can quantize the closed map k into a unitary transformation 
Uh — U K ^ on 7i\ by straightforwardly generalizing the method of [2J. Multiplying this 
unitary operator by the quantum projector 11/ = X/Q-ej \Qj)(Qj\i we obtain the quantum 
open baker's map (j4.38j) 

To obtain an agreement with the notations of £12.4.31 we can interpret the set I = I± U I2 
as the "wall" of the quantum dot, while the complementary interval L = I \ I represents 
the "openings" of the dot, perfectly connected with the "leads". In the previous sections, 
we studied the resonances, that is, the eigenvalues of Bh = U Kt hHi (or of Bh when choosing 
the Walsh quantization). Now, we want to study the "transport" through the dot, using 
the formalism presented in ^2.4.31 We assume that the opening L splits into two disjoint 
"leads" L = L\ U L 2 , and we study the transmission matrix from lead L\ to lead L 2 (for 
simplicity, both leads will have the same width). This matrix is obtained by decomposing 
the scattering matrix (j2.15|) 

sip) = n L J2 ( eM? u >> n ')" <?* u h 

n>0 

into 4 blocks, according to the decomposition Ii^ = © ^-l 2 ■ The transmission matrix 
from L\ to L2 is defined as the block 

(6.1) t(#) = ^ U l 2 U h (Uj U h ) n - 1 U Ll d ^ f ^ e in » tn ■ 

n>l n>l 

Because Tl^ and YIl 2 have the same rank M = N\Li\, £($) is a square matrix of size M. 
According to Landauer's theory of coherent transport, each eigenvalue T»(i?) of the matrix 
£*(■#)£(■#) corresponds to a "transmission channel". The dimensionless conductance of the 
system is then given by the sum over these transmission eigenvalues: 

(6.2) (/(0)=tr(f(tf)t(0)). 

A transmission channel is "classical" if the eigenvalue T« is very close to unity (perfect 
transmission) or close to zero (perfect reflection). The intermediate values correspond to 
the "nonclassical channels" , the importance of which is reflected in the noise power 

(6.3) P0&) = tr (t*tC&) (Id - t m tC&))) , or the Fano factor FC&) = . 

v J g{v) 
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In general it may be necessary to perform the "ensemble averaging" (averaging over i?) to 
obtain significant results [HZ|. However, for the model we will study below, both conduc- 
tance and noise power will depend very little on so this averaging will not be necessary 
To alleviate notations we will supress the dependence in $ in the transmission matrix t. 

Our model. In the remainder of this section, we will compute the quantities characterizing 
the transport through the "dot" when Uh is a Walsh-quantized baker's map similar to the 
operator ()5.9j) . but with D = 4 instead of D = 2. The sequence of values of h consequently 
is given by 

2nh k = A~ k , fc = l,2, 

We will choose the two leads L\ = [0, 1/4] and L 2 = [3/4, 1]: this way, the projectors IL;. 
and III = Id — 11^ — Hl 2 can be represented as tensor product operators: 

= 7r <g> Id 4 (g) • • • ® Id 4 , 
riL 2 = ^3 ® Id 4 ® • • • ® 7d 4 , 
11/ = 7Tf <8> 1^4 ® • • • <2) Id± , 

where 7Tj = |e«) (e^| is a rank-1 orthogonal projector acting on C 4 , and we note 717 = tti ©7r 2 . 

The "inside" propagator for this model, namely = UhUi, is the first one among the 
two quantum maps mentioned in Remark 15.21 its nontrivial spectrum satisfies the fractal 
Weyl law with exponent v = |, and is concentrated near the radius r (B) = 2~ 3 / 4 . 

The number of scattering channels in each lead is the rank of IT^ (equal to that of II^ 2 ). 
It is given by | of the total dimension, and we denote it by 

(6.4) M(h) = ^(2tt/i)- 1 = 4 fc "\ he{h k }. 

The number of channels is "macroscopic" , and each channel is "fully coupled" to the leads. 
We are therefore in a very nonperturbative regime, where resonances have no memory at 
all of the eigenvalues of the closed (unitary) system. 

6.2. Conductance. We will crucially use the fact that all operators under consideration 
act nicely on the tensor product structure Ti^ = (C 4 )®^, that is, they do not entangle the 
quDits. It is then suitable to compute the trace of t*t in a basis adapted to this tensor 
product, and we naturally choose the computational (or position) basis. The conductance 
is then given by 

tr(tt)= (Qj\n\Qj)= E ll^)!! 2 - 
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Let us consider an arbitrary j = eie 2 • • -e^ with e\ = 0, so that < j < 4 fc_1 — 1. Using 
(16. lj) we write 

(6.5) t|g,) = ]T e w n L2 ^ (n, u h y~ l n Ll |q,) = £ e ^ t n \Q,) , 

n>l n>\ 

so that 

ll*IQi)ll 2 = £ e4{n " m) "(^ictnig,). 

m,n>0 

From now on, we replace the notation |Qj), J £ [0,4 fc_1 — 1], by the symbolic notation |e), 
where the sequence e = e 2 ■ • ■ corresponds to j. 

6.2.1. Classical transmission channels. To understand the action of t n on |e) we notice that 
11/ Uh acts on tensor products as 

If n < k, we obtain 

(6.6) t n \e) = n 3 e €n+1 ® e en+2 ® • • • e fc ® ^e ® ^7rje £2 ® • • • ® ^7Tje en , 
frow which we draw the 

Lemma 6.1. Consider a sequence e = 06263 • ■ -e^, and assume that there exists an index 
2 < £ < k such that eg G {0, 3}. Let £ be the smallest such index. Then 



me) 



i/ ^ = 0. 

1 z/ e* = 3 



JTiis shows that |e) is a classical transmission channel. 

Proof. For any 1 < n < £ ~ 2, e n+1 G {1,2} by assumption. Hence the first qu-Dit on 
the right hand side of ()6.6|) vanishes and t n \e) = 0. Furthermore, the state (UiUhY°^ 1 \e) 
admits as first quDit 7T/e Qo = 0, so that t n \e) = for any n > £ . The only remaining 
term in (|6.5|) is t^ _i | e) : 

• if 6e = 0, the first quDit of that term is = 0, so £* _i|e) = t\e) = 0. 

• if e e = 3, ^ -i I e) = e t3 <g> e £fQ+1 g) • • • <g) T*^ ® Fle e2 <%>■•■ T%e H ^ . Since T% is 
unitary, ||^ _i|e)|| = ||t|e)|| = 1. 

□ 

The number of classical channels discussed in Lemma 16.11 is easy to compute: it is 
obtained by removing from the set [0, 4 fc_1 — 1] = { e 2 ■ ■ ■ e k G (Z 4 ) fc_1 } the sequences such 
that eg G {1,2} for all 2 < £ < k (these will be called "nonclassical sequences"). The 
number of the latter is 2 fc_1 , so the number of classical channels is 4 fc_1 — 2 fc_1 . Among 
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them, half are fully reflected, t\e) = 0, and half are fully transmitted, ||t|e)|| = 1. Hence 
the conductance through these classical channels is 

tr d (ft) = . 

Remark 6.1. Such classical channels are mentioned in the analysis of [HZ] for the trans- 
mission through an open kicked rotator. They sit in the phase space regions above the 
lead L\ which are either sent back to L±, or sent to L 2 through the classical dynamics, in 
a time smaller than the Ehrenfest time T^ r = logiV/(log4) = k. For our baker's map B, 
these regions are vertical strips of widths 4 , I = 2, . . . , k which exit to L\ or L 2 at time 
I. The particularity of the Walsh quantization is the exact full transmission (or reflection) 
through these channels. 

6.2.2. Nonclassical transmission channels. The nonclassical channels are necessarily combi- 
nations of the position states |e) with ee G {1, 2} for all 2 < £ < k ("nonclassical" sequences 
or states). The associated positions 4Qj = • e 2 e 3 • • • lie close to the Cantor set C, such 
that r_ = C x I is the set of points never escaping through B or B (see Eq. (13.6(0 . 

For such a state |e), the term ()6.6|) vanishes for all n < k, due to the first qu-Dit 
7i"/e £n+1 = 0. That state therefore accomplishes k "unitary bounces" inside the cavity, 
before it starts to decay out of it. We first consider the terms tk +m \e) for < m < k: 

tk\e) = (e 3 /2) (g) T\e e% ® TX&e 3 • • • ^l&e k , while for < m < k 

t k+m \e) = irzTte em+1 ® ^Je £m+2 ■ • • T\e €h ® Tl{e 3 /2) ® T^iTle^ ■ ■ ■ J^inJ^e^ . 

An explicit computation shows that 

IKa^ejU 2 = ^, lk/-^4 e il| 2 = 2 ' J = <V--,3, 

so that 

(6-8) ||4 +m |e)|| 2 = I ^, 0<m<k-l. 

For larger times n = pk + m, p > 1, m G [0, k — 1], the state t n \e) is obtained from (|6.7|) 
by inserting the operator (vr/jF^) p_1 in front of each quZ)it e er Since irjj-'l has a spectral 
radius |A + | < 1, the norms of these states will decay exponentially fast as n — > oo. This 
argument gives the following 

Lemma 6.2. For any < O < 1, there exists C > such that, for any k > 1 and any 

nonclassical state |e), we have 

m>[efcj 



DISTRIBUTION OF RESONANCES FOR OPEN QUANTUM MAPS 51 

Neglecting errors of order O(2 _0fc//2 ), we just need to compute || Y^m=i^k+ m \^)\\ 2 ■ From 
()6.8|) we already know the diagonal terms: 

Lefcj 

(6-9) ^||t fe+m |e)|| 2 = - + 0(2- efc ). 

m=0 

In the next lemma we will show that the contribution to the conductance of the nondiagonal 
terms is negligible in the semi classical limit. 

Lemma 6.3. Let < O < 1/5. There exists C = C(Q) > such that for any k > 1, 
# { nonclassical e, 3m, m! G [0, Qk], (e\t* k+m t k+m/ \e) ^ } ^ C2 - k / 2 

{ nonclassical e } ~~ 

In other words, in the semi classical limit, a "generic" nonclassical state |e) will satisfy 
Mm, m' G [0, Qk], m ^ m! (e\t* k+m t k+m/ \e) = . 

Proof. Take an arbitrary nonclassical state |e), and any m,m! G [0,OA;], m > m' . From 
(|6.7|) . the first {k — m) quDits of the states t k+m \e) and t k+m i\e)) are respectively 

7T 3 :F 4 *e ero+1 ® J r 4 *e em+2 ® • • • ® ^e,, , 

^3^4e em , +1 <8> ^ 4 *ee m , +2 ® • • • ® -^e, fc+m ,_ m • 

Due to the unitarity of T% and the fact that the ej form an orthonormal basis of C 4 , the two 
states t k+m \e), t k+m i\e) will be orthogonal if the sequences e m+ 2 • • • e k and e m / + 2 • • • ek+ m '-m 
are not equal. Since we took m < Qk, these two sequences are subsequences of length 
(k — m — 1) > (1 — Q)k of the sequence e, shifted from one another by (m — m') steps. 

If the two subsequences are equal, then all subsequences 

e fc-(p+l)A+l • ■ • Cfc-pA ? P = 0, • • • , R , 



A = [m — m ) , H — 



k — m — li 



A 

have to be equal to each other. Hence e contains a subsequence of length (R + 1) A which 
is periodic with period A. 

Let us count the number of such sequences e. Once we have fixed the A bits e k ^A+i • • • e fc, 
the remaining free bits are €i • ■ ■ €k-{R+i)A- The number #(m, m!) of such sequences is 
therefore 2 A x 2 fc_ ( i?+1 ) A_1 . From the definition of R, we get 

#(m,m') < 2 2m - m ' < 2 2ek . 

Taking into account all possible pairs (m, m'), we obtain the following bound for the number 
of nongeneric nonclassical channels: 

# { nonclassical e, 3m,m', < rri < m < Qk, (e\t* k+m t k +m'\e) ^ } < (0A;) 2 2 2efc . 

Since # { nonclassical e } = 2 fc_1 and 2B < 2/5, we have proven the lemma. □ 
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From Lemma l6~21 and Eq. (|6.9jh a generic nonclassical sequence e will satisfy ||i|e) || 2 = ^ + 
(9(2~ 0fe / 2 ). For a nongeneric nonclassical sequence e, we use the simple bound ||t|e)|| 2 < 1. 
As a result, we get the following estimate for the "nonclassical conductance": 

(6.10) tr noncl (t*t) = J2 ll t l e )H 2 + E ll^)f = ^(l + 0(2- efc/2 )). 

nonclassical nonclassical 
generic nongeneric 

Adding this to the "classical conductance" , we get the full conductance 

(6.11) g{&) = tr(f t(0)) = + O(2 (1 - 0/2)fc ) . 

The implied constant is independent of $ G [0, 27r) and < < 1/5. The number of 
scattering channels in our model is given by M(h) = A k ~\ see (IQ1) . so we have proven 
fjl.ljl in Theorem El 

6.3. Noise power. The conductance corresponds to the first moment of the distribution 
of transmission eigenvalues. It can not distinguish between a purely classical transport 
(7* G { 0, 1 }) and a quantum one (some Tj take intermediate values). To do so, we need to 
compute the second moment of these eigenvalues, that is, the trace 

tr((rt) 2 )= Yl rwir, 

or equivalently the noise power ([6.3)1 . As in the previous section, we split the sum on the 
right hand side between the classical and nonclassical states \Qj) = \e). 

Lemma 16.11 shows that half the classical states are in the kernel of t*t, half in the 
eigenspace of t*t associated with the eigenvalue 1 (full transmission). As a consequence, 
the trace over the classical states takes the value 

Ak—l ofc— 1 

trci((t*i) 2 ) = • 

Obviously, the classical states are noiseless. 

The contribution of the nonclassical states is more delicate. According to Lemma 16. 2| 
for any nonclassical state |e) we have (for any < < 1) 

Lefcj 

t\e) = Y,e m H k+m \e) + 0{2~ ek ' 2 ). 

m=0 

We now apply to each state tk+ m \e), m < Qk, the adjoint operator 

(6.12) r = ^e~ ini) r n . 

n>0 

According to ([6.7)1 . the state tk +m \e) has the form 

t k+m \e) = e 3 ® TI-kiWi <S> ■ ■ ■ ® f%niw k , 
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for some explicit set of quDits W£ G C 4 , 2 < £ < k. From the expressions 

Tl I Ul(v l ® ■ ■ ■ <g> v k ) = TiiFiVk ®vx®---® v k ^i , 
we can easily write the action of the operators £* on t& +m |e): 

if n < k, then t* n t k+m \e) = ■K G -K I w k - n+1 ® ... = 0. 
The first non-trivial case of n = k is given by 

t* k t k+m \e) = 7roJ" 4 e3 <S> 7TiW 2 ® - • ■ TrjW* , 
while for any 1 < m' < A; — 1, we have 

7r J r 4 vr/w fc _ m / + i <g> 7T/ J r 4 vr/w;fc-m'+2 <8> • • • niT^iw k ® itiTt&z <g> -niw 2 ® ■■■ niw k - m ' ■ 

The above state is obtained by first evolving |e) A;+m times through the "inside propagator" 
UfJIi, then projecting on the "output lead" L 2 , then evolving backwards (through HiU%) 
k + m' times, and finally projecting on the "input lead" L\. 

As for the case of the operator t, we see that by increasing ml we increase the number 
of quDits on which we apply the operator 7T/JF 4 . Therefore, for any index m, the norm 
of |^w im (e)) will decrease exponentially fast with ml. As in Lemma 16.21 we truncate the 
expansion (|6.12j) to the range m! < Qk, thereby omitting a remainder 0(2~ ek l 2 ). 

We now replace the qu-Dits wg by their explicit values, which depend on the index m. 
We introduce the following notations for operators on C 4 : 

V a p == KaFtftpFl , with a, (3 e {0, J, 3} . 

The explicit decomposition of |Vw,m(e)) depends on the sign of A = f m — m', and on 
whether m,m' = or not (we will often omit to indicate the dependence in e): 

• for m' = m, 

, g 13 . I^o.o) = Voseo ® e t2 ® ■ ■ ■ e ek , 

\ip m ,m) = V ie ® V n e t2 ® ■ ■ • P//e em ® V I3 e tm+1 ® e tm+2 ® ■ ■ ■ e ek . 

• for m = m' + A, A > 0, 

(6.14) 

IV'o.a) = ^03e £A+1 ® e £A+2 ® e efe ® 7rj^Xe ® 717 ^e £2 ® • • • TT/J^e^ 

|Vw,m'+A> = *Poie €A+1 ® V n e €A+2 ® • • ■ V n e €A+mf ® 



54 S. NONNENMACHER AND M. ZWORSKI 

• for m = m' + A, A < 0, 

\^\A\,o) = 7roJ zr 4 e efc __ |A|+1 <g> ^i^e tk _ ]Al+2 ® • • • 71/ J^ 4 e efc <8> P/3e <8> e e2 ® • • • e efc _ |A| , 

(6.15) |VW|A|,m) = 7r 0-^4e efe _| A | +1 ® 7Tj J=ie efc _ |A|+a ® • • • TTjJ~ i&e k ® V n e ® 

® V H e t2 ® • • • <g> P//e em ® V I3 e tm+1 ® e em+2 ® • • • e £fc _ |A| . 

We notice that each of these states contains subfactors e £m+2 <E> • • • <E> e efc if m > m', and 
e em+2 <8> • • • <8> e e fe+m _ m / if m < m'. Compared to its position in the tensor product expansion 
of |e), this subfactor is shifted by m! — m = —A steps. From this remark, and using similar 
methods as for Lemma ESS we can easily prove the following 

Lemma 6.4. Let < 6 < 1/6 and for any pair of indices (m, ml), denote A = m — m! . 
There exists C = C(0) > such that 

#{nonclass. e : 3m 1 ,m' 1 ,m 2 ,m' 2 G [0,6A:], Ai ^ A 2 , ('0 TO ' 1 , mi |V'm 21 m 2 ) ^ °} < £<2~ Cfc 

# { nonclass. e } ~~ 

In other words, for a generic nonclassical sequence e, any two states \if) m ' mi (e)) , |Vw,,m 2 ( e )) 
with mj, m'j < Ofc will be orthogonal to each other if Ai ^ A 2 , that is, if the shifts between 
their respective forward and backward evolution times are different. 

From now on we assume that e is a generic nonclassical sequence in the sense of the 
above lemma. If we group the states |Vv+A,m( e )) into 

|*A(6)) = ^2 l^m+A,m(e)) , 

0<m,m'<tdk 
m=m'+A 

then genericity implies that A(e)|\I / A'( e )) = if A 7^ A'. The square-norm of t*t\e) is 
then given by 

(6.16) ll***|e)U 2 = II^A(e)|| 2 + C(2- efc/2 ). 

|A|<0fe 

As we will see, no further simplification occurs in this expression, meaning that two different 
states |Vw,m) with the same A will generally interfere with each other. Our remaining task 
consists in computing each square norm on the right hand side of (|6.16|) . We will use the 
explicit tensor decompositions (j6.13H6.l5j) . and notice that the overlap between two states 
|Vw,m) is the product of the overlaps of their tensor factors. We split the lengthy, yet 
straightforward computation according to the value of A. 

6.3.1. Norm o/^o- We have 

(6.17) n*oir= Nw*ir +2 M^mW*,*) ■ 

m<0fc 0<m<n<Qk 
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The successive diagonal terms take the values 

||^o,o|| 2 = H^eoH 2 = 777 , while for m > 1 , 

lb 



U m , m \\ 2 = IWojeoinUlWue^) \\V I3 e em+1 \\ 2 = - (jjp ' ' 

1=2 

The sum over the diagonal terms is therefore 



£ ll^ m || 2 = |r + 0((3/8) 0fc ). 
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m<0k 



The nondiagonal terms read, for any < n < Qk, 

n 1 1 n 1 1 

(V^olVVn) = (Vo3e ,Voie ) {X\(e tn Vne tl ))(e en+1 ,V n e en+1 ) = - (-)"" - 



8 y 2' 
1=2 

and forO<m<n<0/c, one similarly gets 

/ , I , > 1 / 3 \ m— 1 1 i ? / 1 \ n— m— 1 1 
\tym,mWn,n) = J (, g J l^' 4 ' 

The sign is + if e m+ i = 1, and — if e m+i = 2. Adding up the real parts of these off-diagonal 
terms, we obtain 

2 J2 M^m,m\^n,n) = ^ + 0(2- @k ). 

0<m<n<ek 

We notice that this contribution is of the same order as the diagonal one. Summing the 
diagonal and nondiagonal parts yields the norm 

21 

(6.18) ll*oir=a„+0(2- efc ). 

oil 

6.3.2. Norm of ^ a with A > 0. From Eq. (j6.14j) we notice that all states \ip m ,m+A), < 
m < Qk — A share the same A last qu-Dits, which results in a common factor 

A 



Y\_ iK/^ej 2 = in the norm ii^aIi 2 - 

1=1 

To avoid taking this factor into account at all steps, we rather consider the states \if}' m m+ A) 
obtained by removing these last A qu-Dits from iVv^m+A)- 
We first compute the square- norm || ip 0;A \\ 2 = A. 

For all m > 0, the first quDit of (VCttm-a) * s ^We £A+1 . From the explicit expression of 
Voi, this qu-Dit vanishes if eA+i = 2, so that 

(6.19) if e A+1 = 2, ||^ A ( e )f = 1-L for all 1< A < Qk . 
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In the opposite case eA+i = 1, the states \ip' mrn+ &) are nontrivial: 



for any < m < Qk — A, 

[Bk]-A 



It ||2 _ £ Am-1 \ 



<V4 m+A l^ „+a) = - t - for 1< m < n < Qk - IAI 



(6.20) so that £ IIC im+ Al| 2 = ^ + 0((3/8) efc - A )) 

m=0 

We then compute the off-diagonal terms: 

— 1 — i 1 1 

(^o, A IC m +A> = T^TT 4 for < m < 0/c - |A|, 

1 Am-1 1±» 1 1 

16 2 n ~ m— ^ 4 

(the sign ± in the last line depends on eA+ m +i)- Summing up the real parts yields: 
2 £ Re(CfAlW = + ^(2- efc+A ) • 

0<m<n<6fc-A 

Adding this to the diagonal terms (jfi.2()j) . restoring the factor 2~ A , and using ()6.19j) results 
in the following norm (which explicitly depends on e): 

(6.21) ||^A(e)ir = ^(^ eA+1=1 + ^ eA+1=2 )+0(2- efc ) forany 1 < A < Qk . 

6.3.3. Norm of ^a with A < 0. As in the previous case, we notice from ()6.15|) that all 
components of \I/a share the same |A| first qu-Dits, which contribute a factor 

a 

( 6 -22) lkoJ r 4e ek _ 1A|+1 || 2 (jJll 7rj - F4ee *-|A|+/ll 2 ) = 4^- 

1=2 

We call Vrra+|A| m ^ ne states with these |A| qu-Dits removed. They have the norms 

Wa|,oII 2 = ^ and formal, ||C + |A|,J| 2 = \ (fr -1 \- 
Hence, the diagonal contribution reads 

L0fcJ-|A| 

E iic + iA,, m ir = ^+0((3/8) efc -i A i). 



m=0 

The nondiagonal terms take the values 



(V>f A |,oK + |A|,J = i J for < n < Qk - |A| 



16 2™" 1 4 

(V>' . iai J< . iai J = - (-)™" 1 — — r - for 1< m < n < Qk - IAI 
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These contributions sum up to 

2 £ Re<C + |A|, m K + |A|,J = ~ + 0(2- efe+ l A l) . 

O<m<n<0fc-|A| 

Putting together this with the diagonal contributions and restoring the factor ()6.22|) . yields 

(6.23) ||vf/ A || 2 = l_L + (2-e*), -Qk < A < — 1 . 

6.3.4. Summing up. We can now sum over all shifts A, |A| < Qk for a given generic 
nonclassical sequence e. The sum over the shifts A < is simple, and independent on the 
sequence e: 

£ ||* A (e)|| 2 = ! + (^2^). 

-e/c<A<o 

The sum over the shifts A > is slightly more delicate, since the norm of |^a(g)) depends 
on e — see Eq. fl6.21|) . However, we notice that the set of generic nonclassical sequences 
can be partitioned into "mirror pairs" (e, e) such that 

We[2,k], e e = l^e e = 2. 

Summing the norms over a "mirror pair" is easy: 

for 1 < A < 0*, ||¥ A (e)|| 2 + ||* A (e)|| 2 = -L i + 0(2^ k ) . 

This contribution is identical (up to the remainder) with ||\EL_ A (e) || 2 + ||^_A(e)|| 2 , which 
shows a sort of symmetry between positive and negative shifts. Summing over all |A| < Qk, 
we get, for any mirror pair (e, e) of generic nonclassical sequences: 

||rt|e)|| 2 + ||rt|e)|| 2 = — + £>(2- efc/2 ) . 

80 

Using Lemma f6.4| we obtain the trace over the nonclassical states: 

tr noncl ((rt) 2 ) = 2 fc - 2 (| + 0(2~ Ck ) + 0{2-W)) . 

Substracting this expression from the "nonclassical conductance" (|6.10|) . and calling C = 
min(C, |0), we finally obtain the noise power: 

(6.24) P = tx(t*t - (t*t) 2 ) = tr noncl (rt - (t*t) 2 ) = 2 k ~ x (H + 0{2~ dk )) . 

This proves ()1.2j) in Theorem El As remarked in the factor 11/80 is close to the 

random- matrix prediction for this quantity, namely 1/8 j^EZl- This is in contrast with 
our remark I3TXI that the semiclassical resonance spectrum of the propagator inside the dot, 
Bh = Ufjli, is quite different from that of a random subunitary matrix. Somehow, the 
matrix t, obtained by summing iterates of B^, has acquired some "genericity" , as far as 
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the distribution of its singular values is concerned. It would be interesting (but quite 
cumbersome) to compute the higher moments of that distribution. 
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